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1.  INTRODUCTION 

In  the  turbomachinery  industry,  the  emphasis  is  to  develop  more  and  more 
efficient  engines  with  higher  reliability  at  a  minimum  cost.  In  the  long  run,  the 
bulk  of  the  cost  is  encountered  when  there  are  compromises  on  reliability  and 
thus  maintainability  of  the  engine.  Millions  and  eventually  billions  of  dollars  can 
be  saved  by  increasing  reliability  and  decreasing  maintenance  down  time  required 
for  the  advanced  turbomachine.  The  primary  cause  for  decreased  reliability  and 
high  maintenance  costs  of  the  engine  is  vibration  during  operation.  The  study  of 
how  to  control  and  minimize  the  vibration  loading  input  to  the  various  members 
of  the  engine  assembly  is  by  no  means  simple 

An  ideal  method  of  solving  these  vibration  problems  is  to  be  able  to  an¬ 
alytically  predict  the  vibrations  during  the  design  and  development  phase  of 
the  engine.  One  tool  used  by  Rotordynamicists  today  is  the  transfer  matrix 
method.  The  method  developed  by  Mykelstad  (1944),  Prohl  (1945),  Koenig 
(1961),  Sevcik  (1963),  Lund  et  al.,  (1967),  Hibner  (1975),  and  Various  other 
authors  is  maintaining  its  popularity  due  to  its  relative  ease  in  use  and  fast 
computer  solution  time. 

This  work  is  mainly  concerned  with  improving  the  accuracy  and  range  of  ap¬ 
plication  of  the  transfer  matrix  method  used  in  rotordynamics  analysis.  Current 
day  rotordynamics  transfer  matrix  programs  are  based  on  modelling  the  rotor 
sections  with  closed  form  beam  theory  solutions.  However,  many  advanced  en¬ 
gines  today  such  as  used  in  aircraft  applications  have  thin  walled  conical  section 

This  thesis  follows  the  style  of  the  ASME  Journal  of  Vibration,  Acoustics, 
Stress,  and  Reliability  in  Design. 
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rotors  which  cannot  be  modeled  as  beam  sections.  In  addition,  the  whole  engine 
must  be  modeled  including  the  housings  and  casings  which  also  do  not  deform 
according  to  beam  theory. 

There  have  been  recent  publications  proposing  various  methods  for  solving 
this  problem  of  inaccuracy  due  to  the  beam  modelling  used  in  the  transfer  ma¬ 
trix  approach.  Rouch  et  al.,  (1979),  Nelson  (1980),  To  (1981),  Greenhill  et  al., 
(1985),  and  Gupta  (1986)  have  all  worked  on  developing  new  finite  elements  for 
the  finite  element  method  as  an  alternative  to  the  transfer  matrix  method.  In  the 
pursuit  of  increased  accuracy,  many  major  turbomachinery  engine  manufactur¬ 
ers  have  generated  highly  complex  and  very  expensive  finite  element  modelling 
of  complete  engine  structures.  Becker  *»♦  ’  .  (1988)  has  published  one  such 

general  finite  element  method  Rr  rotordynamics  analysis.  However  more  accu¬ 
rate  the  analysis  may  become,  there  is  a  high  dollar  cost  for  the  engineering 
and  computer  time  involved  in  the  finite  element  modelling  of  today’s  advanced 
turbomachinerv  such  as  the  light  weight  high  performance  jet  aircraft  engines. 

Shell  finite  element  analysis  was  used  to  generate  the  design  data  presented 
in  this  work.  The  design  data  can  now  be  used  in  the  transfer  matrix  method 
by  design  engineers  and  analysts  without  having  to  do  any  shell  finite  element 
modeling.  Not  only  will  this  result  in  engineering  and  computer  cost  savings 
but  it  will  also  free  up  engineers  and  analysts  to  work  on  other  problems  in 
turbomachinery  which  need  to  be  solved.  Beamish  et  al.,  (1988)  has  presented 
a  paper  on  this  approach,  however,  using  the  finite  element  frequency  method 
of  determining  the  design  data.  The  finite  element  frequency  method  was  also 
published  by  Beamish  et  al.,  (1986).  Thompson  et  al.,  (1988)  modeled  the 
complete  LM2500  U.S.  Navy  gas  turbine  engine  using  the  transfer  matrix  method 
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with  finite  element  conical  section  flexibilities.  Although  the  conical  section 
flexibilities  were  determined  for  only  the  specific  LM2500  gas  turbine  engine, 
his  demonstration  serves  as  an  excellent  validation  of  the  shell  finite  element 
transfer  matrix  approach  presented  in  this  work. 

The  design  data  are  developed  herein  by  modeling  conical  sections  of  in¬ 
dustry  specifications  using  shell  finite  elements.  The  flexibility  coefficients  for 
the  complete  range  of  structural  geometries  encountered  within  industry  have 
been  generated.  These  values  are  tabulated  and  plotted  in  nondimensional  form 
for  ease  in  application  by  turbomachinery  design  engineers  and  analysts.  The 
method  of  replacing  the  old  inaccurate  beam  flexibility  coefficients  with  the  shell 
generated  flexibility  coefficients  adds  no  complexity  to  the  current  transfer  ma¬ 
trix  analysis.  Rather,  it  actually  decreases  computer  solution  time  by  decreasing 
the  number  of  stations  in  the  analysis  and  can  greatly  increase  accuracy  of  crit¬ 
ical  speed  predictions,  especially  in  rotors  with  thin  conical  sections. 

In  this  work,  axisymmetric  cantilevered  shells  with  nonaxisymmetric  loads 
were  modeled.  The  left  end  of  the  structure  was  held  fixed  and  the  right  end 
cantilevered.  The  right  end  was  then  tip  loaded  first  by  a  unit  moment  load 
and  second  by  a  unit  shear  load.  The  tip  radial  and  rotational  displacements 
were  measured  for  the  two  separate  load  cases  in  turn.  By  using  the  force, 
displacement  and  flexibility  mathematics  relations,  the  needed  flexibilities  and 
thus  nondimensional  flexibility  coefficients  were  obtained. 

From  the  beginning  of  this  work,  the  effects  of  cross  section  ovalization  of 
the  hollow  structures  under  load  were  raised  as  a  potential  problem.  The  effects 
on  the  linearity  of  the  flexibility  coefficients  was  found  negligible  for  the  loading 
expected  in  turbomachinery  applications.  A  literature  survey  on  the  theory  of 
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this  ovalization  was  conducted  and  is  presented  including  actual  quantification 
of  these  effects. 
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2.  THEORETICAL  BACKGROUND 

In  order  to  demonstrate  the  improvement  of  using  shell  elements  in  the 
transfer  matrix  method,  a  review  of  the  current  method  is  presented.  The  spe¬ 
cial  considerations  of  using  shell  elements  is  also  presented.  The  field  transfer 
matrix  for  a  general  elastic  structure  is  developed.  The  special  case  of  the  Timo¬ 
shenko  beam  element  using  advanced  strength  of  materials  (A.S.M.)  for  the  field 
transfer  matrix  is  also  presented.  The  flexibilities  of  both  the  general  structure 
element  and  specific  case  beam  element  are  developed.  Convergence  of  flexibility 
values  as  the  general  structure  approaches  a  beam  structure  is  discussed.  The 
flexibility  transformations  of  a  general  structure  loaded  at  opposite  ends  in  turn 
is  developed. 

The  point  transfer  matrix  for  the  conical  shell  structure  cannot  be  modeled 
by  conventional  lumping  at  the  ends  of  discrete  beam  sections.  Large  errors  for 
the  lumped  masses,  polar  and  diametral  moments  of  inertia  will  result  when  the 
structure  is  modeled  in  this  manner.  The  lumped  masses,  polar  and  diametral 
moments  of  inertia  are  shown  as  developed  in  Beamish  et  al.,  (1988). 

The  nonaxisymmetric  partial  differential  equations  for  shells  of  revolution 
are  presented  and  discussed.  The  finite  element  method  of  solution  was  chosen 
because  of  its  ease  in  this  application.  The  shell  finite  elements  are  discussed 
along  with  the  method  of  solution. 
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2.1  Field  Transfer  Matrix  Development 

The  field  transfer  matrix  can  be  developed  for  an  elastic  structure  of  any 
geometry  provided  the  proper  boundary  conditions  are  applied.  The  energy 
and  displacement  principles  of  both  Betti’s  Reciprocal  Theorem  and  Maxwell’s 
Reciprocal  Theorem  (Tuma,  1969)  and  the  theory  of  advanced  strength  of  ma¬ 
terials  have  been  applied  in  the  development  of  the  transfer  matrices  contained 
herein.  Both  ends  of  the  general  elastic  structure  have  been  constrained  to  re¬ 
main  rigid.  These  are  essential  boundary  conditions  in  order  that  the  theory  of 
transfer  matrices  can  accurately  model  the  strain  energy  state  and  deflections 
of  the  structure  under  load.  The  deflections  were  measured  as  shown  in  Fig.  1. 
The  general  form  of  the  transfer  matrix  which  governs  the  relation  between  the 
state  of  the  structure  at  the  left  end,  (i  -  1),  and  at  the  right  end,  (i),  is; 


The  state  variables  of  radial  displacement,  angular  displacement,  moment  load¬ 
ing  and  shear  loading  are  illustrated  in  Fig.  2.  (For  the  derivation  of  equation 
(1),  see  APPENDIX  A.) 
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2.1.1  General  Structure  Element 

It  is  shown  in  APPENDIX  A  for  a  general  structure  element  that  the  state 
variables  at  the  right  end  can  be  related  to  those  at  the  left  by  the  flexibilities  of 
the  element.  This  can  be  stated  in  the  form  of  the  transfer  matrix  of  equation 
(2) 

I  1  l  l  fyO  l(fyo)  ~  fyy  (  ^<-1 

I  I  0  1  fe9  l{fee)~fey  )  #,--i 

\  Mi  l  0  0  1  l  | 

I  1]  )  L0  0  0  1  J  \  l)-i 

It  is  important  at  this  point  to  note  from  equation  (2)  how  the  general  flexi¬ 
bility  coefficients,  a,  6,  a.  and  J,  of  conical  structures  can  be  obtained.  The 
flexibilities  of  the  structure  fye,  fyy,  fee ,  and  f$y,  shown  in  equation  (2),  can  be 
obtained  by  inverting  the  stiffness  matrix  of  the  structure.  Thus  the  problem  of 
generating  the  previously  unknown  conical  shell  structure  transfer  matrix  flex¬ 
ibility  coefficients  is  reduced  to  obtaining  the  inverse  of  the  stiffness  matrix  of 
the  structure.  This  has  been  achieved  by  the  tip  loading  method  explained  in 
the  Introduction.  (For  the  derivation  of  equation  (2),  see  APPENDIX  A.) 
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2.1.2  Timoshenko  Beam  Element 


It  can  be  shown  for  a  Timoshenko  beam  element  that  the  state  variables 
at  the  right  end  can  be  related  to  those  at  the  left  by  the  use  of  closed  form 
expressions.  This  relation  can  be  stated  in  the  form  of  the  transfer  matrix  of 
equation  (3). 
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(3) 


This  closed  formed  transfer  matrix  is  developed  using  advanced  strength  of  ma¬ 
terials.  (For  the  derivation  of  equation  (3),  see  APPENDIX  A.)  Current  day 
transfer  matrix  computational  programs  are  based  on  different  forms  of  equation 
(3).  Some  neglect  the  shear  term  (  qxciC,  )  an(*  others  who  include  the  term  often 
use  a  generalized  constant  shear  form  factor,  K,  =  0.75.  Both  of  these  current 
practices  can  severely  limit  the  applicability  of  transfer  matrix  analysis  to  the 
wide  spectrum  of  rotor  modeling  discussed  in  the  Introduction.  The  errors  are 
quantified  in  section  5.4. 
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2.1.3  Cantilevered  Shaft  Flexibilities 

In  order  to  determine  the  conical  shell  flexibilities  a,  b,  a,  and  3,  it  is  first 
advantageous  to  reduce  these  to  expressions  based  on  the  tip  deflections  of  the 
cantilevered  conical  section  under  unit  tip  loading.  These  are; 


a  =  l(Ym)  +  K, 

(4) 

b  =  Ym 

(5) 

Cl  —  l(Qm)  +  0V 

(6) 

u, 

II 

3 

(7) 

(These  equations  are  developed  in  APPENDIX  B.) 

2.1.4  Cantilevered  Beam  Flexibilities 


The  flexibility  coefficients  a,  6,  a,  and  3  for  a  Timoshenko  beam  element 
including  shear  from  advanced  strength  of  materials  can  be  obtained  using  the 
methods  of  section  2.1.3.  (This  is  also  shown  in  APPENDIX  B.)  The  closed 
form  expressions  are; 


a  = 


b  =  q 
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(9) 

(10) 


2.1.5  Shaft-Beam  Convergence 

The  theory  used  in  the  development  of  equations  (4),  (5),  (6)  and  (i)  has 
no  basic  assumptions  limiting  the  geometry  of  the  structure.  Thus  it  should 
follow  that  as  the  geometry  of  the  structure  approaches  that  of  a  Timoshenko 
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beam  element,  the  corresponding  flexibility  coefficients  a,  b,  a ,  and  J  should 
converge  upon  the  closed  form  expressions  of  equations  (8),  (9)  and  (10).  This  is 
demonstrated  and  discussed  in  the  numerical  results  and  discussion  of  section  5. 

2.1.6  Conical  Flexibility  Transformation 

By  applying  the  principles  of  both  Betti’s  Reciprocal  Theorem  of  energy 
and  Maxwell’s  Reciprocal  Theorem  of  displacements  (Tuma,  1969)  with  transfer 
matrix  mathematics,  a  very  simple  transformation  is  developed  relating  the 
conical  shell  flexibility  coefficients  of  a  positive  angle  cone  to  the  same  geometry 
but  of  a  negative  angle.  (See  Fig.  3.)  The  equations  are; 


—  &p 

(ii) 

bn  = 

(12) 

Q-  n  ~  bp 

(13) 

$n  —  i^p 

(14) 
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(See  APPENDIX  C  for  the  derivation  of  equations  (11)  through  (14).)  For  the 
conical  flexibility  transformation  illustration  of  Fig.  4,  both  the  large  left  end  and 
small  right  end  are  constrained  to  remain  rigid.  The  solid  line  walls  represent  the 
deformed  structure  under  load,  either  shear  or  bending  or  combined.  By  first 
considering  the  undeformed  negative  angle  structure,  fixing  the  left  end,  and 
loading  the  right  end  (positive  attitude  loading),  the  deformed  structure  of  the 
solid  line  results.  Conversely,  by  fixing  the  right  end  of  the  positive  angle  struc¬ 
ture  and  loading  the  left  end  (negative  attitude  loading),  the  identical  deformed 
structure  of  the  solid  line  results.  This  is  due  to  the  rigid  end  constraints.  It  is 
important  to  note  that  the  defformed  structures  of  each  attitude  loading  have 
identical  stress  and  strain  energy  states.  The  tip  deflections  are  measured  as 
previously  defined  in  Fig.  1.  The  equivalent  rotational  displacements  are  easily 
seen.  The  radial  displacement  transformations  are  developed  in  APPENDIX  C. 
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and  negative  attitude 

Unloaded  structure, 
negative  attitude,  (small  end  fixed). 

Unloaded  structure, 
positive  attitude,  (large  end  fixed). 


Fig.  4.  Conical  flexibility  transformation  illustration. 


16 


2.2  Point  Transfer  Matrix  for  Conical  Sections. 

The  point  transfer  matrix  represents  the  lumped  masses,  polar  moments  of 
inertia  and  diametral  moments  of  inertia  of  the  rotor  shaft  or  casing  sections 
lumped  at  the  "point”  between  any  two  sections  (or  stations)  of  the  complete 
rotor-bearing  system.  Other  dynamic  effects  of  the  complete  system  are  assem¬ 
bled  into  the  point  transfer  matrix  but  are  not  discussed  here  because  of  their 
independence  of  the  shaft  modeling.  These  other  dynamic  contributions,  for 
example,  are  bearings,  seals,  supports  and  foundations.  (See  Fig.  5  for  an  illus¬ 
tration  of  the  "point”  transfer  section.)  In  order  to  describe  the  terms  in  the 
"point”  transfer  matrix,  consider  the  simplified  special  case  of  forward  circular 
synchronous  whirl  due  to  unbalance  response; 


1 

0 

0 
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r  Yi-i ) 

0 

1 

0 

0 

j 
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(15) 


When  modeling  a  conical  section  as  a  shell  rather  than  a  series  of  stepped  beams, 
the  lumped  masses,  polar  moments  of  inertia  and  diametral  moments  of  inertia 
cannot  be  evenly  divided  between  the  two  points  on  either  end  of  the  conical 
shaft  section.  Dividing  evenly  is  the  correct  method  for  cylindrical  sections,  but 
large  errors  up  to  changes  in  order  of  magnitude  will  result  if  this  is  attempted 
for  conical  sections,  especially  for  large  cone  angle  sections.  Beamish  et  al., 
(1988),  by  taking  the  geometrical  properties  into  account,  has  developed  closed 
form  expressions  for  the  lumped  masses,  polar  moments  of  inertia  and  diametral 
moments  of  inertia  for  conical  sections. 


17 


R  L 

©  i-1  =  @i-i 


R 


"Point" 

Station 

(1=0) 


"Field" 

Station 

(1=1) 


Fig.  5.  "Point"-"Field"  shaft  sections  diagram. 
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2.2.1  Lumped  Masses  and  Polar  Moments  of  Inertia 

Referring  to  Fig.  6,  according  to  Beamish  et  al.,  (1988),  the  lumped 
at  the  small  radius  cone  end  is 

Cm 


Mi  =  (1 


/ 


-)M 


and  the  lumped  mass  at  the  large  radius  cone  end  is 

\I2  = 

where  the  mass  of  the  entire  conical  section  is 
M  =  P-^-{M0  ~  Min) 

=  ^  i (Rl  -RiR2^  R\)o  -  (R\  -R1R2-  R22)> 


and  the  cone  center  of  mass  is 


Cm  = 


The  outside  radius  is 


Rch0  4 


and  the  inside  radius  is 

Rc 


( Rcm  0)M0-(RCvt  )A  f.n 

'{R2l+2RlR2+3R22y 
{R\  4-  ^^2  +  R\)  ‘ 

{R\^2RxR2^3R\) 


l 


4  [  (R2l+RlR2~R22) 

The  lumped  polar  moment  of  inertia  at  the  small  radius  cone  end  is 

(/„>,  =  (1  -  c-f)ip 

and  the  lumped  polar  moment  of  inertia  at  the  large  radius  cone  end  is 

(iPh  =  (Ci)ip 

where  the  polar  moment  of  inertia  for  the  entire  conical  section  is 


frrd  r 

10 


{R\  +  R*R2  +  R\R\  ^RXR\~  R\) 


{R\  +  R\R2  +  R\R\  +  RXR\  -  R\)in 


mass 

(16) 

(17) 

(18) 

(19) 

(20) 

(21) 

(22) 

(23) 


(24) 


19 


20 


2.2.2  Lumped  Diametral  Moments  of  Inertia 

Referring  to  Fig.  6  and  section  2.2.1,  the  diametral  moment  of  inertia  for 
the  entire  cone  is  given  by 


Id  =  p*' 


I(R\  +  R\R2  -  R\R\  +  R,Rl  +  R\) 


,3,*;  fij. 

30  10  5  ; 
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(R\  -  R\R2  -  R\R\  +  RxR\  -  R\) 


(25) 


3  Rf  R:R2  Rj 

1  30  10  5  ’ 


and  the  lumped  diametral  moment  of  inertia  at  the  small  radius  cone  end  is 
solved  for  by  the  substitution 


l  ~  h  =  Cm 

{Id)\  =  Id  with  ^  (R2)0  =  (RCm)o 

(Rl)tn  =  (Rc„)‘ 


(26) 


and  the  lumped  diametral  moment  of  inertia  at  the  large  radius  cone  end  is 
solved  for  by  substitution 

(  1  ~  l2 

{Id)2  —  Id  with  <  (f?i  )0  =  {Rcxt  )o  (27) 

(  {R\)in  =  (72cm  )in 

The  outer  and  inner  radii  at  the  cone  section  center  of  mass  are 

(flcj,  =  )(  fl2-.fi, )„ 

=  («,).„  ^(~)(Rz-R,hn 


(28) 
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2.3  Shell  Theory 

The  analysis  of  uniform  shells  of  revolution  under  edge  loading  is  by  no 
means  a  trivial  problem.  The  loading  is  nonaxisymmetric,  thus  resulting  in  a 
nonaxisymmetric  distribution  of  stresses  and  strains  throughout  the  shell.  The 
solution  to  the  nonaxisymmetric  partial  differential  equations  for  the  cylinder 
(simplified  special  case  of  cone  with  cone  angle  8  =  0,  see  Fig.  3.)  is  also  nontriv¬ 
ial.  An  exhaustive  literature  search  was  conducted  resulting  in  no  closed  form 
solution  found  for  either  conical  or  cylindrical  shells  of  revolution. 


2.3.1  Shells  of  Revolution 


As  stated  in  Timoshenko  et  al.,  (1959),  the  governing  equations  for  shells 
of  revolution  are 

‘V*.  ,  *0a 


=  -z 

r2 

with  the  nonaxisymmetric  partial  differential  equations 


(29) 


and 


d  ,v  x  dN9,t, 

^(-v*,^)  - — — r j  -  N0jrlcoso,  +  1  rxr3  =  0 


d  c)S'e 

■^■^3^0,0.)  -  ^  y010,ricoso3  t  Xr3r i  =  0 


(30) 


(31) 


(See  Fig.  7  for  an  illustration  of  these  normal  and  shearing  forces.)  It  is  important 
to  note  that  the  solution  of  these  equations  will  yield  the  membrane  stresses  only. 
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2.3.2  Finite  Element  Solution 

The  finite  element  method  was  chosen  to  solve  for  the  flexibility  coefficients 
for  shells  of  revolution.  An  assemblage  of  thin  “quadrilateral  shells”  from  the 
ANSYS  Engineering  Analysis  System  (1985)  were  used.  (See  Fig.  8  for  an  illus¬ 
tration  of  the  ANSYS  quadrilateral  shell  element.)  The  global  stiffness  matrix 
was  assembled  by  ANSYS  taking  into  account  both  membrane  (in-plane)  and 
bending  (out-of-plane)  stresses.  The  element  is  actually  a  thin  flat  plate  derived 
from  thin  plate  theory  simplified  by  the  Kirchhoff  assumptions: 

i)  Straight  lines  initially  perpendicular  to  the  reference  surface  remain 
straight  and  perpendicular  to  that  surface  after  deformation.  (There¬ 
fore.  no  out-of-plane  shear  deformation  is  included  in  the  analysis.) 

ii)  Line  segments  initially  perpendicular  to  the  undeformed  reference  sur¬ 
face  suffer  neither  extension  nor  contraction.  (No  thickness  changes.) 

The  element  stiffness  matrix  is  assembled  by  accounting  for  the  stiffness  contri¬ 
butions  due  to  in-plane  stresses  and  out-of-plane  stresses  separately  and  then 
assembling  these  contributions  together  to  form  the  element  stiffness  matrix. 
First,  for  the  in-plane  membrane  stresses,  from  the  additional  Kirchhoff  assump¬ 
tion: 

iii)  The  transverse  stresses  are  assumed  small  in  comparison  to  in-plane 
membrane  stresses.  (Thus  for  the  strains,  eIZ  =  eyz  =  ezz  =  0). 

and  the  theory  of  elasticity,  the  strain  relations  become; 

€  X  ~ 


du 

dx 

dv 

dy 

du 

Oy 


dv 

dx 


(32) 

(33) 


(34) 
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In  addition,  for  the  in-plane  stresses,  small  rotational  stresses  are  accounted  for 
using  the  method  of  Zienkiewicz  (1977).  Next,  for  the  out-of-plane  bending  finite 
element  shape  functions,  the  above  Kirchhoff  assumption  (iii)  does  not  apply. 
The  method  of  the  DKT  Triangle  of  Batoz  et  al.,  (1980)  and  Razzaque  (1984) 
is  used  for  the  development  of  these  out-of-plane  bending  strain  contributions 
to  the  element  stiffness  matrix. 

Once  the  strain  relations  and  displacement  functions  are  defined  as  above. 
ANSYS  uses  the  energy  method  of  solution  for  the  loaded  structure.  The  prin¬ 
ciple  of  virtual  work  used  is 

SV  =  SU  (35) 

where  the  virtual  strain  energy  is  expressed  as 

6U  =  f  { 6e}T{a}d{Vol )  (36) 

JVol 

where  the  strain  vector  is  as  defined  in  equations  (32),  (33),  (34)  with  the  dis¬ 
placement  functions  of  Zienkiewicz  (1977),  Batoz  et  al.,  (1980)  and  Razzaque 
(1984)  substituted  in  the  equations  for  u  and  v.  The  stress  vector  is  associated 
with  the  strain  vector  by  the  constitutive  equations. 
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3.  OVALIZATION  RESEARCH 

The  ovalization  research  was  accomplished  in  order  to  quantify  the  effect 
of  ovalization  on  the  shell  flexibility  coefficients.  It  was  determined  that  the 
effects  were  negligible  for  the  shell  geometries  and  loading  applicable  in  the 
rotordynamics  industry. 

3.1  Ovalization  of  Shells  of  Revolution 

One  objective  of  this  research  was  to  investigate  and  determine  the  effect 
of  ovalization  on  the  transverse  bending  behavior  of  thin-walled  shells  of  rev¬ 
olution  as  used  for  example  in  gas  turbine  engine  casings.  Ovalization  is  a 
geometric  nonlinear  deformation  process.  This  nonlinear  behavior  to  collapse 
of  shell  structures  under  a  bending  load  is  a  complex  problem  even  for  shells  of 
simple  geometry. 

L.G.  Brazier  (1927)  is  attributed  with  the  first  theoretical  analysis  attempt¬ 
ing  to  quantify  ovalization  in  his  paper  entitled  “On  the  Flexure  of  Thin  Cylin¬ 
drical  Shells  and  other  “Thin”  Sections.”  For  pure  bending,  Brazier  (1927)  used 
the  theory  of  minimum  strain  energy  coupled  with  the  theory  of  flexure  displace¬ 
ment  for  a  circular  cylindrical  thin-walled  shell  of  infinite  length.  To  obtain  the 
equations  of  equilibrium,  first  and  second  order  terms  were  taken  into  account 
giving  rise  to  the  nonlinear  equations.  According  to  Brazier  (1927),  the  shear 
effects  quantified  in  third  and  higher  order  terms  need  not  be  taken  into  account 
as  long  as  the  cylindrical  shells  are  “sufficiently  thin.”  With  the  equations  de¬ 
veloped  for  an  infinitely  long  cylinder,  no  boundary  condition  end  effects  had 
to  be  accounted  for.  Because  maximum  ovalization  occurs  at  maximum  length 
over  radius  ratio  (;)  ,  any  decrease  in  the  (£)  ratio  will  decrease  the  amount 
of  ovalization  predicted  by  Brazier  (1927)  for  an  infinitely  long  cylindrical  shell. 
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As  seen  in  Fig.  9,  the  linear  theory  of  flexure  developed  by  St.  Venant  (as  stated 
in  Brazier  (1927))  predicts  a  much  different  force-displacement  relationship  than 
the  nonlinear  prediction  given  by  Brazier  (1927).  Note  here  that  the  axis  cen¬ 
terline  curvature,  c,  is  inversely  proportional  to  the  radius  of  curvature,  rc.  It 
should  be  noted  that  St.  Venant’s  solution  had  only  taken  first  order  terms  of 
flexure  into  account,  thus  predicting  the  bending  behavior  of  a  non-deforming 
circular  rigid  cross  section  along  the  entire  length  of  the  cylindrical  shell.  How¬ 
ever,  it  should  also  be  noted  that  for  small  displacements  of  cylindrical  shells 
of  small  (j)  ratios,  St.  Venant’s  linear  theory  of  flexure  is  very  accurate.  Theo¬ 
retically,  all  finite  length  cylindrical  shells  will  behave  in  a  manner  somewhere 
between  the  two  limiting  cases  of  St.  Venant  and  Brazier  (1927). 

The  need  here  is  to  quantify  the  effect  of  ovalization  on  the  flexibility  co¬ 
efficients  for  cylindrical  shells  of  practical  engineering  application  such  as  gas- 
turbine  engine  casings.  Referring  to  Fig.  9,  St.  Venant’s  solution  is  represented 
by  the  equation; 

Mb  =  Enr3tc  (37) 


However,  the  governing  differential  equation  presented  by  Brazier  (1927)  is; 


d6i 


d4  v  d2  v 


du >6  di /’4  dil)2 


-18  cV 

~T2 


(1  -  cr2)sin(2xl') 


(38) 


(These  variables  are  illustrated  in  Fig.  10.)  The  solution  is  represented  by  the 
equation; 


Mb  = 


E7rr3t 

2 


2c  - 


3r4c3(l  -  u2) 

? 


For  a  maximum  bending  moment  denoted  by  point  A  of  Fig.  9; 

—  _  2  v'2£,7rrf2 
~~  9(1  -</2)? 


(39) 


(40) 


and, 
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Au>  =  (41) 

(See  Fig.  11  for  an  illustration  of  this  Aio  Brazier  flattening  (ovalization).)  It 
should  be  noted  that  for  pure  bending,  the  Brazier  flattening  of  equation  (41) 
gives  the  maximum  ovalizing  at  point  A  of  Fig.  9.) 

According  to  Antonenko  (1981),  '‘'‘The  phenomenon  of  ovalization  of  a  cir¬ 
cular  cross  section  in  the  bending  of  a  shell  as  a  beam  is  attributable  to  the  effect 
of  additional  loads  which  arise  in  the  process  of  curvature  of  the  shell  axis.  The 
pattern  of  loading  may  be  described  as  follows.  In  bending  of  the  shell  axis,  ten¬ 
sile  axial  forces  on  the  convex  side  and  compressive  forces  on  the  concave  side 
yield  resultants  directed  toward  the  shell  axis.  Ovalization  of  the  cross  section 
takes  place  under  the  influence  of  these  loads.”  (Refer  to  Fig.  12  for  a  graphical 
representation  of  these  resultant  forces.)  Antonenko  (1981)  accounted  for  oval¬ 
ization  in  an  infinite  system  of  coupled  differential  equations  which  he  solved 
by  successive  approximations.  Antonenko  claims  that  ovalization  need  be  taken 
into  account  only  for  cylindrical  shells  with  large  (£)  ratios  and  very  small  (£) 
ratios  (thin-walled)  only.  No  limiting  ratios  were  given.  However,  limiting  ratios 
are  established  in  the  following  Ovalization  Research  Section. 


Curvature,  C 
(0  increasing) 


Fig.  9.  Brazier's  nonlinear  solution  compared  with 
St.  Venant's  linear  solution. 


Fig.  1U.  Brazier  s  coordnate  system  for  a  cylindrical  cross 
section. 
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C:  Compressive  force 
T:  Tension  force 
R:  Resultant  force 
D.E.:  Differential  element 


Fig.  12.  Antonenko  mechanization  of  ovalization. 
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3.2  Computational  Analysis 

Stephens  et  al.,  (1975)  generated  finite  difference  models  using  the  STAGS 
computer  code.  The  models  were  isotropic  cylindrical  shells  with  (£)  =  .01, 
E  =  1.0  x  I07psi.  and  u  =  0.3.  The  (£)  ratio  was  varied  from  model  to  model. 
A  moment  couple  for  pure  bending,  A/,  was  applied.  Both  ends  of  the  cylindrical 
shell  were  essentially  rigid  (no  ovalization  at  ends).  The  finite  difference  grid 
point  spacing  was  varied  in  both  the  axial  and  circumferential  directions  to 
allow  for  a  finer  mesh  concentration  at  the  expected  maximum  ovalization  region 
(midpsan  compression).  The  STAGS  code  used  a  linear  flexure  theory  coupled 
with  a  nonlinear  geometric  “fudge  factor’’  to  allow  for  the  expected  nonlinear 
prebuckling  state. 

The  deformation  shape  output  by  the  above  mentioned  STAGS  computer 
code  was  controlled  to  deform  as  shown  in  Fig.  11.  Fig.  13  illustrates  the  oval¬ 
ization  dependence  upon  the  (£)  ratio.  As  shown  in  equation  (42),  the  critical 
bending  moment,  A/cr,  is  solely  dependent  upon  the  cylindrical  shell  geometry 
and  Young’s  Modulus  of  elasticity. 

Mcr  =  O.GETrrt2  (42) 


M  is  the  bending  moment  applied  given  in  equations  (39)  and  (40).  Stephens 

et  al.,  (1975)  combined  equations  (40)  and  (42)  to  obtain  the  moment  ratio  for 

a  Brazier  infinite  cylindrical  shell  for  this  particular  shell  model  geometry; 

M  _  2y2£'7rrf2  ^  1 

MCr  9(1 -//2)  5  0  .GEnrt2 

or  for  a  poisson’s  ratio  =  0.3,  by  substitution  and  algebra; 


(44) 
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Also,  by  combining  the  model  geometry  of  £  =  .01  and  Brazier’s  equation  (41), 

Atu  2  r  ^ 

-  =  -  x  -  =  22.22  (45 

t  9  t 

It  is  important  to  note  in  Fig.  13  that  the  envelope  curve  for  the  model  ap¬ 
proaches  the  predicted  value  of  Brazier’s  equation  for  an  infinite  cylindrical  shell. 
It  is  also  important  to  note  that  for  all  (£)  ratios,  ovalization  does  occur  but 
is  negligible  for  small  (£)  ratios.  Circumferential  flattening  or  ovalizing  of  the 
cross  section  is  more  and  more  nonlinear  for  higher  and  higher  (~)  ratios.  The 
major  percentage  of  ovalization  (when  large  ovalization  does  occur)  takes  place 
at  or  near  the  buckling  load.  Refer  to  Table  1  for  calculations  of  actual  change 
in  diameter  over  thickness,  (^)  ratios,  for  the  models  represented  in  Fig.  13. 
It  is  very  important  to  note  from  Table  1  that  for  small  bending  moments  and 
neutral  axis  displacements,  ovalization  becomes  very  small. 

Additional  computational  analysis  was  performed  for  cylindrical  and  conical 
sections  using  ANSYS.  The  results  confirm  negligible  ovalization  for  geometries 
and  loadings  of  practical  interest.  (See  the  numerical  results  and  discussion  of 
section  5.) 


Fig.  13.  STAGS  computer  code  envelope  curve  for  Bra¬ 
zier  flattening. 
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Table  1.  Nondimensional  Ovalization  for  Thin  Cylindrical  Shells  Under 

Varying  Bending  Load  with  (£)  =  .01,  Varying  (£)  ratio. 
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@  Buckling 
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@  \  Buckling 
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3.3  Experimental  Verification 

Experiments  were  carried  out  by  Brazier  (1927)  with  long  (j  =  45)  thin 
tubular  beams  made  of  celluloid  and  loaded  in  pure  bending.  All  specimens 
deformed  into  the  quasi-oval  form  of  Fig.  11  at  high  loads  near  buckling.  The 
nature  of  the  force  verse  neutral  axis  curvature  is  shown  in  Fig.  14.  It  is  very  im¬ 
portant  to  note  that  deviation  from  the  linear  St.  Venant’s  solution  is  negligible 
until  near  buckling  loads  are  applied.  Notice  as  stated  earlier,  that  the  exper¬ 
imental  tubes  do  behave  in  between  the  two  limiting  cases  of  St.  Tenant  linear 
theory  and  Brazier's  (1927)  geometric  nonlinearity  theory.  These  experiments 
by  Brazier  (1927)  also  agree  with  the  computational  results  of  Stephens  et  al., 
(1975). 

Experiments  by  Lundquist  (1933  and  1935)  show  that  for  hundreds  of  metal 
tubes  ranging  with  (£)  ratios  up  to  5  and  very  thin- walled  ratios  (~)  down  to 
.0015,  there  was  no  measured  ovalization.  St.  Venant’s  linear  theory  of  flexure 
was  perfectly  suitable  for  these  specimens.  Osgood’s  (1938)  experiments  record 
no  ovalization  for  metal  specimens  with  (£)  ratios  up  to  10,  but  for  (~)  ratios 
greater  than  10  ovalization  was  observed  however  small.  Roark’s  (1941)  exper¬ 
iments  record  no  ovalization  for  metal  specimens  with  (£)  ratios  ranging  up  to 
2.  Wilson's  et  al.,  (1941)  experiments  record  no  ovalization  for  shells  with  (£) 
ratios  ranging  up  to  2.4.  Gerard  et  al.,  (1968)  concluded  on  these  experiments 
that  the  Brazier  effect  is  only  seen  in  cylindrical  shells  with  high  (-)  ratios, 
thin-walled,  and  even  then  only  at  high  loads  near  buckling. 

According  to  Bushnell  (1981),  the  nonlinear  prebuckling  effects  of  Brazier 
flattening  have  a  significant  impact  upon  the  buckling  load.  However,  at  loads 
well  below  the  buckling  load,  the  effects  are  insignificant.  The  nonlinear  pre- 
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buckling  effects  of  Brazier  flattening  do  cause  a  drop  off  in  cylindrical  shell  stiff¬ 
ness.  However,  this  nonlinear  stiffness  effect  is  only  significant  at  high  Brazier 
flattening. 


Curvature.  C 
(0  increasing) 


Fig.  14.  Brazier's  experimental  plot  compared  with  Bra¬ 
zier  s  theoretical  plot  compared  with  St.Venant  s  linear 
plot. 
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3.4  Ovalization  Conclusions 

The  following  conclusions  are  emphasized  as  the  important  ideas  which 
should  be  understood  from  this  ovalization  study. 

A  closed  form  linear  solution  relating  bending  moment  to  curvature  of 
neutral  axis  does  exist  for  cylindrical  shells  in  pure  bending.  This  is  called 
St.  Tenant's  theory  of  flexure.  This  linear  solution  assumes  no  change  in  cross 
section  throughout  deflection. 

A  closed  form  nonlinear  solution  relating  bending  moment  to  curvature  of 
neutral  axis  does  exist  for  cylindrical  shells  in  pure  bending.  This  is  called 
Brazier’s  theory  of  flexure  which  accounts  for  second  order  terms  introduced  by 
“Brazier  flattening”  or  ovalization. 

The  Brazier  flattening  is  more  pronounced  on  the  compression  side  of  the 
cylindrical  shell. 

Thin- walled  cylindrical  shells  with  low  (^)  ratios  bend  according  to  the 
linear  solution  on  which  the  flexibility  calculations  of  this  work  is  based. 

In  applications  such  as  gas  turbine  engine  casings,  and  rotor  shaft  sections, 
stiff  boundary  conditions  usually  restrict  cross  section  ovalization  thus  making 
the  contribution  of  the  Brazier  flattening  even  more  negligible. 

It  is  very  important  to  note  that  deviation  from  the  linear  St.  Tenant’s 
solution  is  negligible  until  near  buckling  loads  are  applied.  This  is  true  for  even 
the  large  (£)  ratios  and  small  (£)  ratios  of  less  than  .01. 

Because  of  this  Brazier  flattening  deformation  behavior,  for  small  loads  and 
consequently  small  neutral  axis  displacements,  ovalization  of  the  cylindrical  shell 
cross  section  has  almost  no  effect  upon  the  linearity  of  the  load-displacement 
(stiffness)  curves  for  even  large  (£)  ratios  in  transverse  bending. 
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3.5  Additional  Comments 

It  should  also  be  noted  that  from  the  Stephens  et  ah,  (1975)  computational 
solutions,  ambient  pressure  has  an  effect  on  the  Brazier  flattening.  External  pres¬ 
sures  applied  to  a  cylindrical  shell  accentuate  or  increase  the  Brazier  flattening 
while  internal  pressures  lessen  the  Brazier  flattening.  Thus  for  turbomachinery 
with  high  internal  pressure  such  as  in  a  gas  turbine  engine,  even  less  Brazier 
flattening  will  occur.  Brazier  flattening  is  also  affected  by  the  Young’s  modulus 
of  the  material.  In  the  limiting  case,  a  completely  strain  hardened  material  will 
have  a  much  smaller  neutral-axis  deflection  with  load  and  thus  a  much  smaller 
Brazier  flattening. 

The  Brazier  flattening  phenomenon  will  have  an  effect  on  the  dynamical 
behavior  of  cylindrical  shells.  However,  this  effect  will  be  negligible  in  cylindrical 
shells  with  boundary  conditions  limiting  the  Brazier  flattening.  Because  the 
Brazier  flattening  deformation  is  opposite  in  direction  to  the  dynamic  response 
deflection,  a  small  reduction  (if  any)  in  the  first  bending  natural  frequency  of  the 
cylindrical  shell  would  be  expected.  With  the  cyclic  deformation  shapes  expected 
similar  to  Fig.  11,  the  small  fluctuations  in  the  effective  inertia  should  occur  but 
are  expected  to  be  negligible.  Also,  with  Brazier  flattening,  this  increase  in 
deflection  of  the  shell  should  increase  the  differential  element  velocities  and  thus 
increase  the  internal  material  damping,  an  effect  expected  to  be  negligible. 
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4.  FINITE  ELEMENT  COMPUTATIONAL  ANALYSIS 

The  finite  element  method  was  used  to  solve  for  the  flexibility  coefficients 
of  the  conical  structures  of  interest.  The  ANSYS  engineering  analysis  computer 
code  was  used  to  generate  the  shell  models,  boundary  conditions,  meshes  and 
solutions  for  these  conical  structures.  The  models  and  boundary  conditions  are 
defined.  The  meshes  for  these  models  are  discussed  and  sample  illustrations 
given. 

4.1  Models  Defined 

Two  hundred  and  seventy  five  finite  element  models  were  generated.  These 
included  all  combinations  of  the  parameters; 

l 

-  =  .25.  .5,  1.  2,  4 

r 

-  =  .0125.  .025,  .05,  .1,  .2 

r 

9  =  -75°.  -60°,  -45°,  -30°,  -15°.  0°.  15°,  30°,  45°.  60°,  75° 

Each  of  the  two  hundred  seventy  five  had  to  be  tip  loaded,  first  with  a  unit 
shear  load  and  solved,  and  then  with  a  unit  pure  bending  moment  and  solved. 
Five  hundred  and  fifty  solutions  were  generated.  (See  the  sample  schematic  of 
varying  (-)  ratio  and  varying  cone  angle  9  in  Fig.  15.) 
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4.2  Boundary  Conditions 

For  an  illustration  of  the  cantilevered  conical  section  with  loading  applied, 
see  Fig.  16.  It  is  important  to  note  that  all  models  were  fixed  at  the  left  end.  The 
right  end  was  constrained  to  remain  infinitely  rigid  by  restricting  the  degrees  of 
freedom  for  these  end  nodes.  Because  the  free  end  was  constrained  rigid,  the  load 
acted  as  if  it  was  uniformly  distributed  to  each  node  thus  ensuring  no  distorted 
loading  or  local  deformations.  This  is  essential  in  order  to  obtain  the  proper 
displacements  and  thus  flexibility  coefficients.  (See  APPENDIX  D  for  sample 
input  data  for  ANSYS  displacement  constraints  and  rigid  region  constraints.) 

4.3  Mesh  Analysis 

All  models  were  generated  using  the  ANSYS  mesh  generation  procedure. 
In  general,  it  is  important  to  keep  the  individual  elements  of  the  mesh  within 
certain  dimensional  constraints.  We  can  define  these  constraints  as  the  aspect 
ratio,  (^-)  and  the  arc  angle,  0C.  (Refer  to  Fig.  17  for  an  illustration  of  these 
parameters.)  These  constraints  are  important  for  any  finite  element  simply 
because  of  the  limitations  on  the  finite  element  shape  functions  derived  from  the 
particular  structural  theory.  In  order  to  determine  the  specific  limitations  on 
the  aspect  ratio  and  arc  angle,  meshes  of  the  extreme  geometry  models  (( -)  =  4 
with  cone  angle  =  0°)  and  ((£)  =  0.25  with  cone  angle  =  75°)  were  meshed 
first  varying  only  arc  angle  from  30°  to  15°  and  then  varying  only  the  aspect 
ratio  from  15  to  1.  The  solutions  converged  to  within  0.1%  of  the  coarser  mesh 
solution.  The  limits  on  the  mesh  aspect  ratios  were  thus  established  at  0C  = 
15^OJ.  and  (<=#-)  =  6maj;.  (See  Fig.  18  through  Fig.  27  for  samples  of  meshes 
used  for  the  solution  generatin'  See  APPENDIX  D  for  sample  ANSYS  input 
for  model  geometry  and  mesh  ration.) 
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Fig.  18.  ANSYS  mesh  with  (  )  =  1  and  9  =  (J1 
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1  (L./R )  = 


Fig.  21.  ANSYS  mesh  with  (  ‘-)  =  1  and  9  -  45' 
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Fig.  22.  ANSYS  mesh  with  (£•)  =  1  and  B  —  60'  . 


Fig.  23.  ANSYS  mesh  with  (  )  =  l  and  0 


Fig.  24.  ANSYS  mesh  with  (;)  =  '25  and  6  =  45'. 
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Fig.  25.  ANSYS  mesh  with  I '-)  =  .5  and  0  =  45'. 


Fig.  26.  ANSYS  mesh  with  ( '- )  =  2  and  6  —  45 
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Fig.  27.  ANSYS  mesh  with  =  4  and  a  =  45°. 
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5.  NUMERICAL  RESULTS  AND  DISCUSSION 

Ovalization  has  been  demonstrated  to  be  negligible  using  ANSYS  for  trans¬ 
verse  loading  conditions  found  in  turbomachinery.  A  discussion  of  these  oval¬ 
ization  tests  is  presented.  The  flexibility  coefficients  are  presented  in  nondimen- 
sionalized  tabular  and  plot  form  for  ease  in  turbomachinery  design  and  analysis 
application.  The  shell  finite  element  models  have  been  verified  by  compari¬ 
son  with  solutions  from  advanced  strength  of  materials.  The  satisfaction  of  the 
mathematical  flexibility  transformation  identities  is  also  presented  as  verification 
of  the  model  boundary  conditions.  A  discussion  of  current  methods  of  including 
shear  in  the  closed  form  beam  theory  solutions  is  presented. 

Rotor  transfer  matrix  demonstrations  were  performed  on  the  60°  cone  test 
piece  presented  hy  Beamish  et  al.,  (1988)  and  the  Space  Shuttle  Main  Engine 
High  Pressure  Oxidizer  Turbopump  (SSME  HPOTP)  as  presented  by  Li  (1979). 
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5.1  Negligible  Ovalization 

Based  on  the  ovalization  study,  three  models  most  likely  to  ovalize  were 
modeled  and  analyzed.  The  parameters  chosen  for  study  were  (£)  =  4,  (1)  = 
0.0125,  and  cone  angles  0°.  75°.  and  —75°.  A  maximum  load  based  on  the 
worst  case  vibrational  displacement  of  200  mils  (.2  inches)  was  used.  A  load 
vector  based  on  a  linear  extrapolation  was  used  to  check  if  the  nonlinear  anal¬ 
ysis  would  give  the  same  results  as  the  linear  analysis.  The  load  vector  was 
obtained  from  multiplying  the  unit  load  vector  by  200  mils  and  dividing  by  the 
tip  displacement  at  unit  load.  The  large  displacement  nonlinear  solution  option 
within  ANSYS  was  selected  and  each  solution  converged  within  three  iterations. 
The  tip  and  midlength  ( |)  nodal  displacements  were  output.  For  all  three  mod¬ 
els,  there  were  no  changes  in  midlength  diameters  and  the  tip  deflections  were 
200  mils  indicating  a  linear  flexibility  throughout  this  loading  range.  Because 
no  ovalization  was  measured,  all  275  models  were  solved  with  a  linear  analysis 
with  one  load  step 
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5.2  Design  Data  Generation 

The  data  generated  for  the  models  described  in  the  previous  section  were 
nondirnensionalized  for  ease  in  turbomachinery  design  and  analysis  application. 
Tables  2  through  4  list  the  nondimensional  flexibility  coefficients  for  cylindrical 
sections.  The  values  are  listed  for  varying  (4)  and  ( ^ )  ratios.  Tables  5  through 
24  also  list  the  nondimensional  flexibility  coefficients  for  conical  sections.  The 
vaiues  are  listed  for  varying  cone  angle  0,  and  varying  (  - )  and  (-)  ratios.  Fig¬ 
ures  28  through  47  are  plots  of  the  nondimensional  flexibilities  listed  in  Tables  5 
through  24  respectively.  The  plots  may  be  used  for  quick  interpolation  of  nondi- 
rnensional  flexibilities.  Xote  for  the  (aEr^g  =  aEr_o)  and  (a£V2_^  =  bEr2_0) 
tables  and  plots,  that  there  are  negative  values.  The  curves  represent  the  flexi¬ 
bility  characteristics  dominated  by  the  shear  loading.  However,  for  each  curve, 
the  slope  is  changing  and  eventually  as  ( p )  increases  the  flexibility  coefficients 
will  be  dominated  bv  the  bending  moment  loading  and  thus  become  positive. 
Likewise,  the  flexibility  coefficients  (bEr2_0  -  aEr2_$)  and  ( 3Er =  3Er3_e) 
will  always  be  positive  because  each  is  uncoupled  and  represents  the  flexibility 
characteristics  of  the  structure  solely  due  to  pure  bending.  In  pure  bending, 
there  is  no  she;  r  taken  into  account.  It  is  important  to  note  for  all  these  nondi¬ 
mensional  flexibility  plots,  that  the  curves  only  represent  a  best  fit  between  the 
computational  data  points.  The  data  points  were  marked  on  the  plots  with 
symbols  which  represent  the  only  true  computational  data,  the  curves  between 
are  not  computational  data.  The  Tables  of  the  nondimensional  flexibilities  are 
contained  herein  to  give  the  designer  or  analyst  the  option  of  using  a  numerical 
interpolation  scheme  of  his  choice. 
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Table  2.  Nondimensional  Flexibilities  ( aEr )  for  Cylindrical 
Sections  (0  =  0)  with  Varying  (£)  and  (£)  Ratios. 


1 

r 

f  r 

l 

r 

aEr 

0.25 

-16.3875  j 

i 

0.5 

-32.5050 

!  0.0125 

1.0 

-61.8909 

1 

| 

2.0 

-98.0523 

4.0 

. 

10.8891 

0.25 

-8.01212 

0.5 

-16.1382 

|  0.025 

1.0 

-30.7997 

2.0 

-48.8386 

4.0 

5.64177 

r 

0.25 

-3.68059 

0.5 

-7.88378 

0.05 

1.0 

-15.2192 

2.0 

-24.1826 

_ 

4.0 

3.12114  j 

0.25 

-1.39222 

0.5 

-3.61497 

0.1 

1.0 

-7.38548 

2.0 

-11.8026 

_ -  -  - _ 

4.0 

1.94488 

0.25 

-0.354660 

1 

j 

0.5 

-1.35794 

0.2 

1.0 

-3.35416 

1 

• 

2.0 

-5.55745 

i 

Table  3.  Nondimensional  Flexibilities  ( bEr 2  =  aEr 2)  for 
Cylindrical  Sections  (0  =  0)  with  Varying  £  and  -  Ratios. 


t 

r 

1 

r 

bEr 2  =  aEr 2 

1 

0.25 

0.753797 

0.5 

3.13300 

0.0125 

1.0 

12.7363 

| 

2.0 

51.3101 

4.0 

205.622 

0.25 

0.370897 

0.5 

1.54250 

0.025 

1.0 

6.32531 

2.0 

25.5801 

4.0 

102.742 

0.25 

0.184290 

0.5 

0.754806 

0.05 

1.0 

3.13337 

2.0 

12.7302 

_ 

4.0 

51.2805 

0.25 

0.0919208 

0.5 

0.371620 

0.1 

1.0 

1.54363 

2.0 

6.32142 

4.0 

25.5507 

0.25 

0.0458154 

0.5 

0.184268 

0.2 

1.0 

0.756213 

2.0 

3.12746 

4.0 

12.6917 

Table  4.  Nondimensional  Flexibilities  (0Er3)  for  Cylindrical 
Sections  (0  =  0)  with  Varying  j  and  £  Ratios. 


1 

i  t 

r 

- — - 

I 

r 

- ! 

0Erz 

0.25 

6.03040 

1 

0.5 

12.5320 

|  0.0125 

1.0 

25.4725 

1 

| 

2.0 

51.3100 

! - 

4.0 

102.811 

0.25 

2.96719 

i 

0.5 

6.16999 

0.025 

1.0 

12.6506 

2.0 

25.5801 

4.0 

51.3709 

r 

i 

0.25 

1.47433 

1 

l 

0.5 

3.01922 

;  o.o5 

1.0 

6.26673  ! 

2.0 

12.7301 

1 

4.0 

25.6404 

i 

0.25 

0.735368 

0.5 

1.48648 

0.1 

1.0 

3.08726 

2.0 

6.32142 

4.0 

12.7754 

0.25 

0.366524 

0.5 

0.737069 

0.2 

1.0 

1.51242 

2.0  | 

3.12745 

4.0 

6.34587 

Table  5.  Nondimensional  Flexibilities  ( aEr )  for  Conical  Sec¬ 
tions  with  ((£)  =  .0125  ),  Varying  0,  and  Varying  (£)  Ratio. 


0 

l 

r 

aEr+&  —  aEr_& 

0.25 

-16.4088 

0.5 

-31.6478 

15 

1.0 

-57.8755 

2.0 

-93.3205 

4.0 

-104.709 

0.25 

-17.6177 

0.5 

-33.1231 

30 

1.0 

-58.8090 

2.0 

-95.4185 

4.0 

-135.138 

r 

0.25 

-20.5644 

0.5 

-37.6571 

45 

1.0 

-65.2536 

2.0 

-106.495 

4.0 

-165.270 

0.25 

-27.1032 

0.5 

-48.1376 

60 

1.0 

-81.7094 

2.0 

-134.742 

4.0 

-221.627 

0.25 

-45.4020 

0.5 

-77.6723 

75 

1.0 

-130.452 

2.0 

-221.181 

4.0 

-385.649 
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Table  6.  Nondimensional  Flexibilities  ( aEr )  for  Conical  Sec¬ 
tions  with  ((£)  =  .025  ),  Varying  0,  and  Varying  (£)  Ratio. 


0 

l 

r 

aEr+®  =  aEr-@ 

i 

0.25 

-8.02446 

0.5 

-15.7143 

15 

1.0 

-28.8052  j 

2.0 

-46.4997 

4.0 

-52.1826 

0.25 

-8.62022 

0.5 

-16.4428 

30 

1.0 

-29.2559 

2.0 

-47.5162 

4.0 

-67.3361 

0.25 

-10.0659 

0.5 

-18.6749 

45 

1.0 

-32.4185 

2.0 

-52.9506  | 

4.0 

-82.2487 

0.25 

-13.2575 

0.5 

-23.7970 

60 

1.0 

-40.4472 

2.0 

-66.7359 

4.0 

-109.971 

0.25 

-22.0413 

0.5 

-37.8917 

75 

1.0 

-63.6207 

2.0 

-107.963 

4.0 

-189.506 

Table  7.  Nondimensional  Flexibilities  ( aEr )  for  Conical  Sec¬ 
tions  with  ((£)  =  .05  ),  Varying  0,  and  Varying  (£)  Ratio. 


0 

i 

r 

aEr+e  =  aEr_& 

0.25 

-3.68986 

0.5 

-7.67802 

15 

1.0 

-14.2391 

2.0 

-23.0452 

4.0 

-25.8297 

0.25 

-3.97299 

0.5 

-8.03363 

30 

1.0 

-14.4451 

2.0 

-23.5077 

4.0 

-33.3074 

0.25 

-4.64925 

0.5 

-9.11334 

45 

1.0 

-15.9531 

2.0 

-26.0858 

4.0 

-40.5182 

0.25 

-6.12105 

0.5 

-11.5506 

60 

1.0 

-19.7273 

2.0 

-32.5355 

4.0 

-53.6437 

0.25 

-10.0532 

0.5 

-17.9039 

75 

1.0 

-29.9758 

2.0 

-50.6176 

4.0 

-89.4359 

65 


Table  8.  Nondimensional  Flexibilities  ( aEr )  for  Conical  Sec¬ 
tions  with  ((£)  =  .1  ),  Varying  0,  and  Varying  (~)  Ratio. 


0 

l- 

r 

aEr  ,  ~  =  aEr  _ 

tW  — 

0.25 

-1.40197 

0.5 

-3.52347 

15 

1.0 

-6.91415 

2.0 

-11.2706 

4.0 

-12.5724 

0.25 

-1.52682 

0.5 

-3.69121 

30 

1.0 

-6.99747 

2.0 

-11.4472 

4.0 

-16.1763 

0.25 

-1.81070 

0.5 

-4.18548 

45 

1.0 

-7.67514 

2.0 

-12.5752 

4.0 

-19.4612 

0.25 

-2.40248 

0.5 

-5.26514 

60 

1.0 

-9.31562 

2.0 

-15.3136 

4.0 

-25.0916 

0.25 

-3.89256 

0.5 

-7.83074 

75 

1.0 

-13.2275 

2.0 

-21.9353 

4.0 

-38.3099 

Table  9.  Nondimensional  Flexibilities  ( aEr )  for  Conical  Sec¬ 
tions  with  ((£)  =  .2  ),  Varying  0,  and  Varying  (£)  Ratio. 


0 

l 

r 

aEr+e  =  aEr~® 

0.25 

-0.365955 

0.5 

-1.32964 

15 

1.0 

-3.14373 

2.0 

-5.33993 

4.0 

-5.87906 

0.25 

-0.422952 

1 

0.5 

-1.40656 

30 

1.0 

-3.17270 

2.0 

-5.37712 

4.0 

-7.53184 

0.25 

-0.537534 

0.5 

-1.60990 

|  45 

1.0 

-3.44683 

2.0 

-5.78792 

4.0 

-8.84485 

0.2! 

-0.755183 

0.5 

-2.02251 

60 

1.0 

-4.06754 

2.0 

-6.73290 

4.0 

-10.7921 

0.25 

-1.26172 

0.5 

-2.89676 

75 

1.0 

-5.25712 

2.0 

-8.43630 

4.0 

-13.9820 

i 

Table  13.  Nondimensional  Flexibilities  (bEr^@  =  aEr^& 
Conical  Sections  with  ((£)  =  .!,  Varying  0,  and  Varying  ( 


Ratio. 


0 

l 

r 

bEr^@  =  aEr^& 

0.25 

0.407040 

0.5 

1.26040 

15 

1.0 

3.21266 

2.0 

8.05920 

4.0 

19.9812 

0.25 

0.851053 

0.5 

2.45556 

30 

1.0 

5.52912 

2.0 

11.8138 

4.0 

24.9129 

0.25 

1.69223 

0.5 

4.61720 

45 

1.0 

9.77256 

2.0 

19.5604 

4.0 

39.0475 

0.25 

3.86121 

0.5 

9.98283 

60 

1.0 

20.3417 

2.0 

39.7651 

78.6287 

0.25 

13.5442 

0.5 

32.7980 

75 

1.0 

65.4812 

2.0 

128.894 

4.0 

262.191 
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Table  14.  Nondimensional  Flexibilities  (bEr^_@  =  aEr^&)  for 
Conical  Sections  with  ((£)  =  .2  ,  Varying  0,  and  Varying  (£) 
Ratio. 


— 

0 

l 

r 

i 

bErle  =  &Er*@ 

0.25 

0.0820812 

0.5 

0.471078 

15 

1.0 

1.46323  ! 

2.0 

3.85164 

4.0 

9.7002. 

0.25 

0.156430 

0.5 

0.871381 

30 

1.0 

2.44760 

2.0 

5.49738 

4.0 

11.7580 

0.25 

0.333677 

0.5 

1.60694 

45 

1.0 

4.20855 

1 

2.0 

8.80371 

i 

4.0 

17.7554 

0.25 

0.824413 

0.5 

3.38390 

60 

1.0 

8.33288 

2.0 

16.8576 

| 

4.0 

33.5125 

0.25 

2.92445 

0.5 

10.1182 

75 

1.0 

23.3119 

2.0 

46.7298 

4.0 

94.4523 

L 


Table  15.  Non  dimensional  Flexibilities  (aEr^_&  =  bEr^&)  for 
Conical  Sections  with  ((£)  =  .0125  ,  Varying  0,  and  Varying 
(£)  Ratio. 


0 

—  1  ' 

1 

r 

aErl&  =  bEr^_& 

0.25 

-3.35914 

0.5 

-4.81764 

15 

1.0 

-4.03825 

2.0 

4.52193 

4.0 

26.9966 

i 

0.25 

-7.95148 

0.5 

-12.0540 

30 

1.0 

-14.7028 

2.0 

-12.7436 

4.0 

-5.68461 

0.25 

-15.0169 

0.5 

-21.5003 

45 

1.0 

-25.5206 

2.0 

-25.0119 

4.0 

-21.3568 

0.25 

-29.8036 

0.5 

-38.6405 

60 

1.0 

-42.4161 

2.0 

-41.2734 

4.0 

-38.2425 

0.25 

-77.3563 

0.5 

-85.9213 

75 

1.0 

-85.6796 

2.0 

-82.0337 

4.0 

-78.6086 
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Table  16.  Nondimensional  Flexibilities  (a£V^0  =  bEr2_&)  for 
Conical  Sections  with  ((£)  =  .025  ,  Varying  0,  and  Varying 
(~)  Ratio. 


l 

r 

aErl@  =  bEr2_@ 

0.25 

-1.62549 

0.5 

-2.38807 

15 

1.0 

-1.99987 

2.0 

2.27569 

4.0 

13.5125 

0.25 

-3.85713 

0.5 

-5.96493 

30 

1.0 

-7.29022 

2.0 

-6.31742 

4.0 

-2.80081 

0.25 

-7.29076 

0.5 

-10.6205 

45 

1.0 

-12.6356 

2.0 

-12.3952 

4.0 

-10.5965 

0.25 

-14.4526 

0.5 

-19.0011 

60 

1.0 

-20.9050 

2.0 

-20.3661 

4.0 

-18.9281 

-37.0794 


Table  17.  Nondimensional  Flexibilities  (aEr^Q  =  bEr^_@)  for 
Conical  Sections  with  ((£)  =  .05  ,  Varying  0,  and  Varying 
Ra 


0.707560 

1.14735 

0.963623 

1.17096 

6.78615 


1.70794 

2.87817 

3.56147 

3.08032 

1.33143 


3.25015 

5.11992 

6.15827 

6.04743 

5.16653 


-6.44910 

-9.09183 

10.0791 

-9.82485 

-9.15410 


16.2538 
19.1525 
19.1269 
18.3326 
17.9714 
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Table  18.  Nondimensional  Flexibilities  =  bEr^Q)  for 

Conical  Sections  with  ((£)  =  .1  ,  Varying  0,  and  Varying  (£) 
Ratio. 


I  0 

l 

r 

aErl&  =  bEr2_@ 

! 

0.25 

-0.207751 

| 

0.5 

-0.479667 

15 

1.0 

-0.430284 

2.0 

0.635722 

4.0 

3.44023 

0.25 

-0.551703 

0.5 

-1.25374 

30 

1.0 

-1.67476 

2.0 

-1.44016 

| — 

4.0 

-0.571204 

0.25 

-1.09435 

0.5 

-2.24837 

45 

1.0 

-2.88920 

2.0 

-2.84164 

4.0 

-2.40827 

0.25 

-2.21721 

0.5 

-3.96207 

60 

1.0 

-4.62488 

2.0 

-4.50137 

4.0 

-4.17636 

0.25 

-5.47422 

0.5 

-7.86734 

75 

1.0 

-8.03047 

2.0 

-7.59648 

4.0 

-7.44044 

Table  20.  Nondimensional  Flexibilities  (/ 3Er 3)  for  Conical 
Sections  with  ((£)  =  .0125  ,  Varying  0.  and  Varying  ^  Ratio. 


0.25  |  7.14700 

0.5  |  13.5022 


15  1.0  23.2778 

2.0  35.7056 

4.0  47.5998 


30 

0.25 

0.5 

1.0 

2.0 

4.0 

12.6024 

21.7229 

32.9972 

43.5092 

50.4334 

0.25 

28.3401 

0.5 

44.4925 

45 

1.0 

60.5462 

2.0 

72.0681 

4.0 

77.9095 

0.25 

84.1246 

0.5 

118.426 

60 

1.0 

144.894 

2.0 

159.603 

4.0 

165.637 

0.25 

459.849 

0.5 

557.759 

75 

1.0 

611.539 

2.0 

635.678 

_ 

4.0 

644.932 

78 


Table  21.  Nondimensional  Flexibilities  (/ 3Er 3)  for  Conical 
Sections  with  ((^)  =  .025  ,  Varying  0,  and  Varying  £  Ratio. 


f3Erl&  =  /3ES@ 

3.50333 

6.64343 

11.5468 

17.7759 

23.7601 

6.13577 

10.6836 

16.3384 

21.6171 

25.1346 


13.7237 

21.8517 

29.9033 

35.7143 

38.7508 


40.5056 

57.9004 

71.1896 

78.7122 

82.1097 


218.002 

267.920 

295.104 

308.752 

316.701 
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Table  22.  Nondimensional  Flexibilities  (flEr3)  for  Conical 
Sections  with  ((|)  =  .05  ,  Varying  0,  and  Varying  £  Ratio. 


0 

1 

r 

/3Er3+Q  =  0Er3_@ 

0.25 

1.70392 

0.5 

3.23745 

15 

1.0 

5.70020 

2.0 

8.81801 

4.0 

11.8266 

0.25 

2.87902 

0.5 

5.16903 

30 

1.0 

8.01938 

2.0 

10.6646 

4.0 

12.4495 

0.25 

6.27776 

0.5 

10.4957 

45 

1.0 

14.5735 

2.0 

17.4931 

i 

4.0 

19.0701 

0.25 

18.1984 

0.5 

27.4899 

60 

1.0 

34.2474 

2.0 

38.0483 

4.0 

39.9561 

0.25 

94.8717 

0.5 

122.469 

75 

1.0 

136.188 

2.0 

143.171 

1  . 

4.0 

149.201 

Hi 


Table  23.  Nondimensional  Flexibilities  ( 0Er 3)  for  Conical 
Sections  with  ((^)  =  .1  ,  Varying  0,  and  Varying  £  Ratio. 


0.797174 

1.56152 

2.78250 

4.34747 

5.85535 


1.19742 

2.40369 

3.85445 

5.18680 

6.08543 


2.39152 

4.73766 

6.88336 

8.35940 

9.15980 


6.57586 

12.0412 

15.7164 

17.6314 

18.6126 


32.2798 

49.8608 

57.4500 

60.6482 

63.6868 
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Table  24.  Non  dimensional  Flexibilities  (/ 3Er 3)  for  Conical 
Sections  with  ((£)  =  .2  ,  Varying  0,  and  Varying  -  Ratio. 


0.418580 

0.991108 

1.76059 

2.44602 

2.89589 


0.662584 

1.77247 

3.00468 

3.78917 

4.18643 


1.59028 
4.18545 
6.46271 
7.50258 
7.95741 


7.14131 

15.4369 

20.4806 

22.0550 

22.9947 


Fig. 28.  Nondimensional  flexibility  plots  (aEr)  for  coni¬ 
cal  sections  with  ( ‘7)  =  0.0125,  varying  9  and  varying  (  '-) 
tio. 


3 


2 

i 


Fig.  29.  N'ondimensional  flexibility  plots  (a£>)  for  conical 
sections  with  ( t)  =  0.025,  varying  0  and  varying  ( ratio. 


HtH 


Fig.  37.  Nondimensional  flexibility  plots  ( bEr =  a 
for  conical  sections  with  =  0.2.  varying  9  and  varying 


Fig.  39.  Nondimensional  flexibility  plots  ( aEr =  bEr\t) 
for  conical  sections  with  ( l-)  =  0.025,  varying  0  and  varying 


ig-  40.  .Nondimensional  flexibility 


plots  (a£rlH 
varying  9  and 


Fig.  42.  Nondimensional  flexibility  plots  (a£V^  =  6fr2j 
for  conical  sections  with  =  0.2,  varying  9  and  varying 
( ~ )  ratio. 


Fig.  46.  Nondimensional  flexibility  plots  (JFr3)  for  con¬ 
ical  sections  with  ( t)  =  o.l,  varying  B  and  varying  ( ; ) 
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5.3  Model  Verification 

With  any  finite  element  computational  analysis,  the  model  must  be  veri¬ 
fied  in  some  manner  to  ensure  acceptable  computational  results.  The  two  most 
common  methods  are  to  try  and  reproduce  experimental  data  or  try  and  repro¬ 
duce  some  theoretical  closed  form  solution  data.  The  finite  element  models  used 
herein  were  first  verified  by  comparison  to  results  of  theoretical  closed  form  so¬ 
lutions.  The  models  were  then  verified  with  experimental  data  on  the  60°  cone 
test  piece  described  by  Beamish  et  al.,  (1988).  This  is  discussed  in  the  rotor 
transfer  matrix  demonstration  of  section  5.5.1.  For  additional  verification,  the 
mathematical  conical  flexibility  transformations  of  equations  (11)  through  (14) 
were  fully  satisfied. 

5.3.1  Cantilevered  Cylinders 

Six  structures  of  cylindrical  cross  section,  constant  midwall  radius  (r  =  1.0 
inch),  Young’s  Modulus  of  30  x  106  psi,  shear  modulus  of  11.54  x  106  psi,  and 
with  varying  (£)  and  (£)  ratios  as  listed  in  Table  25  were  modeled  with  the 
previously  described  finite  element  procedure.  The  radial  (Y,,)  and  rotational 
(dv)  deflections  for  a  unit  shear  load  were  calculated  by  the  closed  form  ad¬ 
vanced  strength  of  materials  beam  theory  equations  found  in  APPENDIX  A. 
These  values  were  tabulated  for  comparison  with  ANSYS  generated  solutions  in 
Table  25.  One  important  preliminary  concern  was  to  observe  how  the  ANSYS 
finite  element  models  would  behave  for  the  relatively  thick  shells  of  (£)  =  0.2. 
Normally,  thin  plates  or  shells  based  on  Kirchhoff  adhoc  assumptions  (described 
previously)  have  the  limitation  of  (^)  >  .05.  However,  as  listed  in  Table  25,  for 
( ; )  =  20  the  percentage  differences  for  (>'’,.)  and  (0t.)  are  only  1.412  and  1.583 
respectively.  For  these  structures,  the  A. S.M.  solutions  are  expected  to  be  close 
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Table  25.  Cantilevered  Cylinder 

A. S.M.  method  with  F.E.M. 

Displacements 

Comparison 

I 

Beam  Theory 

Shell  Theory 

1 

r 

t 

r 

A. S.M. 

S. F.E.M. 

Difference 

{Y,) 

(Y,) 

(%) 

20 

.0125 

2.3049  E-3 

2.3209  E-3 

.689 

10 

.0125 

3.0363  E-4 

3.0670  E-4 

1.001 

1 

.0125 

;  2.3530  E-6 

2.4876  E-6 

5.411 

20 

.2 

1.4247  E-4 

1.4451  E-4 

1.412 

10 

.2 

|  1.8710  E-5 

1.9035  E-5 

1.707 

1 

.2 

1.3764  E-7 

1.3701  E-7 

-.460 

T 


l 

r 

i 

!  t 

r 

Beam  Theory 

A. S.M. 

m 

, 

Shell  Theory 

S. F.E.M. 

(0v) 

Difference 

(%)  1 

20 

1  .0125 

1.6976  E-4 

1.7139  E-4 

1 

.951 

10 

.0125 

4.2440  E-5 

4.2848  E-5 

.952 

1 

.0125 

4.2440  E-7 

4.2454  E-7 

.033 

20 

.2 

1.0505  E-5 

1.0674  E-5 

1.583 

10 

.2 

2.6263  E-6 

2.6651  E-6 

1.456 

1 

.2 

2.6263  E-8 

2.5207  E-8 

-4.189 
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to  exact.  Thus,  the  thick  shell  finite  element  model  was  accurate.  However,  it 
must  be  noted  that  the  out-of-plane  additional  flexibilities  admitted  by  these 
finite  elements  (discussed  previously)  allow  even  the  thick  shells  to  be  modeled 
accurately.  As  expected,  the  thin  shells  for  (£)  =  20  had  even  better  results 
of  0.689%  and  0.951%  differences  of  the  ( V,. )  and  (0V)  displacements  respec¬ 
tively.  It  should  also  be  noted  at  this  time  that  for  very  small  (£)  ratios,  neither 
A.S.M.  closed  form  solutions  nor  shell  finite  elements  will  be  accurate  due  to  the 
theory  geometrical  limitations.  For  these  geometries,  three  dimensional  solid 
isoparametric  brick  elements  based  on  elasticity  theory  should  be  used.  The 
limits  on  the  very  small  (£)  ratios  have  not  been  established.  However,  for  the 
(£)  and  (£)  ratios  used  in  this  work  with  (/)  always  greater  than  (f),  the  shell 
elements  with  in-plane  and  out-of-plane  considerations  are  quite  accurate. 

Beam  theory  closed  form  solutions  were  generated  for  the  complete  set 
of  model  geometries  for  cylindrical  sections.  The  nondimensional  flexibilities 
( bEr 2  =  aEr 2)  are  listed  in  Table  26.  The  nondimensional  flexibilities  ( 3Er 3) 
are  listed  in  Table  27.  The  beam  theory  closed  form  solutions  are  plotted  to¬ 
gether  with  the  S.F.E.M.  solutions  from  tables  3  and  4  on  figures  48  through 
57.  For  pure  bending,  the  beam  theory  closed  form  solutions  are  expected  to  be 
close  to  exact.  It  can  be  seen  by  the  comparison  of  these  plots  that  the  ANSYS 
models  are  very  accurate. 
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Table  26.  Beam  Theory  Closed  Form  Nondimensional  Flexi¬ 
bilities  ( bEr 2  =  aEr2)  for  Cylindrical  Sections  (©  =  0)  with 
Varying  (£)  and  (£)  Ratios. 


t 

r 

1 

r 

bEr2  =  aEr2 

0.25 

0.795744 

0.5 

3.18298 

0.0125 

1.0 

12.7319 

2.0 

50.9276 

4.0 

203.711 

0.25 

0.397826 

0.5 

1.59130 

0.025 

1.0 

6.36521 

2.0 

25.4608 

4.0 

101.843 

0.25 

0.198820 

0.5 

0.795278 

0.05 

1.0 

3.18111 

2.0 

12.7245 

4.0 

50.8978 

0.25 

0.0992239 

0.5 

0.396895 

0.1 

1.0 

1.58758 

2.0 

6.35033 

i 

4.0 

25.4013 

0.25 

0.0492435 

0.5 

0.196974 

0.2 

1.0 

0.787896 

1 

1 

2.0 

3.15159 

! 

4.0 

12.6063 

I 


Fig.  48.  Beam  theory  closed  form  nondimensional  flex¬ 
ibility  plots  compared  to  finite  element  nondimensional 
flexibility  plots  { bEr 2  =  a£Y3)  for  cylindrical  sections  (9  = 
0)  with  varying  ( l- )  ratios  and  (;)  =  .0125. 


bEr 2  =  aEr 2 
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Fig.  49.  Beam  theory  closed  form  nondimensional  flex¬ 
ibility  plots  compared  to  finite  element  nondimensional 
flexibility  plots  (bEr2  =  aEr2)  for  cylindrical  sections  (9  = 
0)  with  varying  ( \ )  ratios  and  (;)  =  025. 
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Fig.  50.  Beam  theory  closed  form  nondimensional  flex¬ 
ibility  plots  compared  to  finite  element  nondimensional 
flexibility  plots  ( bEr 2  =  aEr2)  for  cylindrical  sections  (0  = 
0)  with  varying  (;)  ratios  and  (;)  =  05. 


bEr3  —  aEr 
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Fig.  51.  Beam  theory  closed  form  nondimensional  flex¬ 
ibility  plots  compared  to  finite  element  nondimensional 
flexibility  plots  (bEr2  =  aEr2)  for  cylindrical  sections  (9  = 
0)  with  varying  (;)  ratios  and  (£•)  =  .1. 
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Fig.  53.  Beam  theory  closed  form  nondimensional  flex¬ 
ibility  plots  compared  to  finite  element  nondimensional 
flexibility  plots  (. 3Er 3)  for  cylindrical  sections  (6  =  0)  with 
varying  ( ; )  ratios  and  ( f )  =  .0125. 


llE  r3 
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Fig.  54.  Beam  theory  closed  form  nondimensional  flex¬ 
ibility  plots  compared  to  finite  element  nondimensional 
flexibility  plots  (,J£V3)  for  cylindrical  sections  (0  =  0)  with 
varying  (  )  ratios  and  ( )  =  .025. 
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Fig.  55.  Beam  theory  closed  form  nondimensional  flex¬ 
ibility  plots  compared  to  finite  element  nondimensional 
flexibility  plots  (JEr3)  for  cylindrical  sections  (9  =  0)  with 
varying  ( l- )  ratios  and  ( \ )  =  .05. 
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Fig.  57.  Beam  theory  closed  form  nondimensional  flex¬ 
ibility  plots  compared  to  finite  element  nondimensional 
flexibility  plots  ( 3Er 3)  for  cylindrical  sections  {6  -  0)  with 
varying  (£)  ratios  and  (~)  =  .2. 


117 


5.3.2  Conical  Flexibility  Transformations 
An  additional  method  of  model  verification  is  the  satisfaction  of  the  math¬ 
ematical  flexibility  transformations  of  APPENDIX  C, 

an  =  ap-  l(bp )  4-  l2(j3p )  -  l(ap) 
bn  =  l{3p )  -  bp 

(CTO) 

ocn  —  /(t^p)  cip 

^ p 

or  more  specifically  for  an  axisymmetric  conical  structure, 


<Zn  —  dp 

(11) 

bn  dip 

(12) 

Otn  —  bp 

(13) 

3n  =  3p 

(14) 

These  mathematical  expressions  would  not  be  satisfied  if  there  were  any  bound¬ 
ary  condition  errors  in  the  models.  This  is  explained  by  the  energy  theorems  of 
Betti  and  Maxwell  (Tuma  (1969)).  All  of  the  flexibilities  of  the  275  models  were 
input  into  the  above  equations  of  (CIO)  and  (11)  through  (14).  All  mathematical 
identities  were  satisfied. 
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5.4  Including  Shear  in  Beam  Theory 

Different  practices  accounting  for  shear  are  currently  being  used  in  rotor- 
dynamic  transfer  matrix  programs.  Some  do  not  include  any  shear  effects  at 
all,  and  most  use  a  constant  shear  form  factor,  K,  =  0.75.  Neither  of  these 
practices  is  adequate  over  the  full  range  of  section  geometries  and  can  result 
in  large  percentage  errors  in  the  flexibility  coefficient  ( aEr ).  Table  28  lists 
the  nondimensional  flexibilities  (aEr)  of  cylindrical  sections  for  the  closed  form 
beam  theory  solutions  using  the  constant  shear  form  factor.  Table  29  lists  the 
nondimensional  flexibilities  (aEr)  of  cylindrical  sections  for  the  closed  form  so¬ 
lutions  with  the  improvements  of  advanced  strength  of  materials  (A.S.M.).  (Ks 
is  determined  by  equation  (A9)).  Table  2  (discussed  previously)  lists  the  nondi¬ 
mensional  flexibilities  (aEr)  of  cylindrical  sections  determined  by  the  shell  finite 
element  analysis.  Figures  58  through  62  illustrate  the  large  percentage  errors 
of  assuming  a  constant  K,  =  0.75.  It  can  be  seen  from  these  plots  that  the 
assumption  of  a  constant  shear  form  factor  can  cause  very  significant  errors  in 
flexibilities.  When  modelling  cylindrical  beam  sections  for  rotor  dynamics  analy¬ 
sis.  The  A.S.M.  closed  form  solution  should  always  be  used.  However,  the  closed 
form  solution  using  A.S.M. is  also  not  exact.  Whether  the  shell  flexibilities  are 
closer  to  the  real  values  than  the  A.S.M.  method  is  unproven  for  these  cylindrical 
sections.  However,  it  is  assumed  that  the  shell  flexibilities  for  the  thin  sections 
are  more  accurate  than  the  A.S.M.  flexibilities.  Because  the  nature  of  the  curves 
don't  change  appreciably  with  thickness,  it  is  further  assumed  that  the  shells 
are  more  accurate  over  the  entire  range  of  section  geometries  contained  in  this 
work. 


Table  29.  Beam  Theory  Closed  Form  Non  dimensioned 
Flexibilities  ( aEr )  for  Cylindrical  Sections  (0  =  0)  with 
Varying  (£)  and  (£)  Ratios. 


Fig.  58.  Beam  theory  closed  form  nondimensional  flex¬ 
ibility  plots  compared  to  finite  element  nondimensional 
flexibility  plots  ( aEr )  for  cylindrical  sections  (0  =  0)  with 
varying  ( ; )  ratios  and  (7)  =  .0125. 


Fig.  59.  Beam  theory  closed  form  nondimensional  flex¬ 
ibility  plots  compared  to  finite  element  nondimensional 
flexibility  plots  (aEr)  for  cylindrical  sections  (0  =  0)  with 
varying  ( -)  ratios  and  (M  =  .025. 
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Fig.  60.  Beam  theory  closed  form  nondimensional  flex¬ 
ibility  plots  compared  to  finite  element  nondimensional 
flexibility  plots  (aEr)  for  cylindrical  sections  (©  =  0)  with 
varying  ( ; )  ratios  and  ( ; )  =  .05. 
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Fig.  61.  Beam  theory  closed  form  nondimensional  flex¬ 
ibility  plots  compared  to  finite  element  nondimensional 
flexibility  plots  (a£Y)  for  cylindrical  sections  (0  =  0)  with 
varying  (  l- )  ratios  and  ( '  )  =  .1. 
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Fig.  62.  Beam  theory  dosed  form  nondimensional  flex¬ 
ibility  plots  compared  to  finite  element  nondimensional 
flexibility  plots  ( aEr )  for  cylindrical  sections  (0  =  0)  with 
varying  ( £ )  ratios  and  (;)  =  .2. 
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5.5  Rotor  Transfer  Matrix  Demonstrations 

The  impact  of  modeling  conical  sections  with  the  shell  flexibility  coefficients 
in  place  of  the  stepped  discrete  beam  flexibility  coefficients  is  significant.  Two 
models  were  picked  to  demonstrate  the  impact  of  using  shell  elements  in  place 
of  beam  elements.  The  first  modeled  was  the  60°  cone  test  piece  of  Beamish 
et  al.,  (1988).  The  fundamental  free-free  natural  frequency  was  determined  and 
the  mode  shape  plotted. 

The  second  transfer  matrix  demonstration  was  performed  on  the  Space 
Shuttle  Main  Engine  High  Pressure  Oxidizer  Turbopump  (SSME  HPOTP).  The 
model  station  data  was  as  reported  in  Li  (1979).  The  SSME  HPOTP  was 
modeled  first  using  beam  element  and  second  using  the  new  shell  element  data 
contained  in  this  thesis  work.  Campbell  diagrams  were  then  plotted  for  both 
runs.  The  critical  frequencies  of  the  two  methods  were  then  compared.  A 
representative  mode  shape  was  also  plotted  for  comparison. 
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5.5.1  60°  Cone  Test  Piece 

The  60°  cone  test  piece  of  Beamish  et  al.,  (1988)  is  illustrated  in  Fig.  63.  The 
experimental  test  fundamental  mode  natural  frequency  was  788  hertz.  Beamish 
et  al.,  (1988)  reported  3231  hertz  when  modeled  with  infinite  stepped  beams  and 
791  hertz  when  modeled  with  shell  flexibilities  by  the  frequency  method.  For  this 
report,  the  cone  was  first  modeled  with  three  stepped  discrete  A.S.M.beams  as 
shown  in  Fig.  64.  The  station  data  was  input  to  the  Kelly  (1989)  transfer  matrix 
rotordynamic  analysis  program.  The  fundamental  first  mode  natural  frequency 
by  this  stepped  beam  method  was  2378  hertz.  The  cone  was  then  modeled 
with  the  new  shell  flexibility  coefficients.  The  shell  flexibility  coefficients  were 
obtained  by  interpolating  directly  off  the  plots  of  design  data  figures  28  through 
47.  The  station  data  was  again  input  to  the  Kelly  (1989)  program  but  with  the 
new  shell  flexibility  coefficients.  The  fundamental  natural  frequency  by  this  shell 
element  transfer  matrix  method  was  734  hertz.  The  corresponding  fundamental 
mode  shapes  for  both  methods  along  with  the  test  mode  shape  of  Beamish  et 
al.,  (1988)  are  plotted  in  Fig.  65.  There  are  large  differences  between  the  natural 
frequencies  predicted  using  discrete  beams  and  the  natural  frequencies  predicted 
using  the  shell  element  data  of  this  report.  This  shows  that  the  beam  modeling 
is  inadequate  for  conical  shell  sections. 


Fig.  63.  60°cone  test  piece  with  dimensioning 
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Fig.  65.  Cone  test  piece  fundamental  mode  shape.  (Test 
frequency  =  788  Hz.  Stepped  beams  frequency  =  2379 
Hz.  Shell  element  transfer  matrix  method  frequency  = 
734  Hz. ) 
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5.5.2  Space  Shuttle  Main  Engine  High  Pressure 
Oxidizer  Turbopump  (SSME  HPOTP) 

The  SSME  HPOTP  illustrated  in  Fig.  66  was  modeled  using  the  station 
data  as  reported  in  Li  (1979).  It  was  first  modeled  as  A.S.M.  beam  elements 
input  to  the  Kelly  (1989)  program.  The  model  was  run  for  eight  different  rotor 
running  speeds  as  shown  by  the  dat?  points  of  the  Campbell  diagram  of  Fig.  67. 
The  SSME  HPOTP  was  then  modeled  with  the  cylindrical  shell  flexibilities 
interpolated  off  the  plots  of  figures  48  through  62.  The  model  was  run  for 
the  same  eight  rotor  running  speeds  as  the  previous  beam  elements  model.  A 
Campbell  diagram  was  plotted  for  these  runs  as  shown  in  Fig.  68.  In  Table  30, 
the  critical  speeds  are  listed  as  predicted  by  both  models.  It  is  important  to 
note,  for  the  low  natural  frequencies,  no  significant  difference  between  the  two 
methods.  However,  it  is  very  important  to  note  that  as  shown  in  sections  5.3 
and  5.4,  there  are  no  significant  differences  between  the  A.S.M.  flexibilities  and 
S.F.E.M.  flexibilities  for  cylindrical  sections.  Both  models  here  are  based  on  the 
Li  (1979)  station  data.  Li  (1979)  had  reduced  the  conical  sections  to  cylindrical 
sections.  If  the  SSME  HPOTP  was  modeled  with  conical  sections  as  conical 
sections  with  the  shell  flexibilities,  differences  in  predicted  critical  speeds  would 
be  expected.  However,  the  impact  on  the  predicted  critical  speed  for  this  specific 
engine  using  S.F.E.M.  in  place  of  the  beams  has  not  been  determined.  The  mode 
shapes  of  the  9,286  rpm  critical  speed  are  plotted  for  both  the  A.S.M.  and  the 
shell  finite  element  transfer  matrix  methods  on  Fig.  69.  Again,  for  the  cylindrical 
sections  and  lower  critical  speeds,  the  A.S.M.  method  is  quite  accurate. 
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Fig.  66.  Space  Shuttle  Mam  Engine  High  Pressure  Oxi¬ 
dizer  Turbopump  Assembly  (SSME  HPOTP). 


RUNNING  SPEED  (rpm) 


Fig.  67. 


Campbell  diagram  for  the  A.S.M.  method. 


RUNNING  SPEED  (rpm) 


Fig.  68.  Campbell  diagram  for  the  shell  element  transfer 
matrix  method. 


Table  30.  SSME  HPOTP  Critical  Speeds  Comparison 
between  the  A. S.M.  Method  and  the  Shell  Element 
Transfer  Matrix  Method. 
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Critical 

Speed 

Critical 

Beam  Theory 

A. S.M. 

Critical 

Speed 

(rpm) 

Shell  Theory 
S.E.T.M. 

Critical 

Speed 

(rpm) 

1 

9,286 

9,286 

2 

10,000 

10,000 

3 

1 

11,116 

10,870 

4 

1  * 

11,785 

11,786 
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6.  CONCLUSION 

Design  data  in  tabular  and  plot  form  were  developed  using  shell  finite  ele¬ 
ments.  The  design  data  are  in  nondimensionalized  form  for  ready  and  easy  use 
by  turbomachinery  design  engineers  and  analysts.  The  ovalization  of  thin  and 
thick-walled  cylindrical  and  conical  engine  sections  was  determined  to  have  a 
negligible  effect  on  rotordvnamics  applications.  The  nondimensionalized  design 
data  generated  are  functions  of  conical  section  axial  length,  small  end  radius, 
wall  thickness,  cone  angle  and  Young’s  modulus  of  elasticity.  The  nondimension¬ 
alized  parameters  of  these  sections  cover  the  complete  range  of  rotors,  casings, 
and  housings  used  throughout  the  turbomachinery  industry.  The  design  data 
are  presented  in  the  form  of  flexibility  coefficients  which  should  be  input  to  the 
transfer  matrix  analysis.  The  need  for  these  flexibility  coefficients  to  be  used  in 
the  current  day  transfer  matrix  computer  programs  has  been  demonstrated. 

It  is  highly  recommended  that  these  design  data  be  incorporated  into  the 
existing  computer  codes  for  transfer  matrix  analysis.  The  computer  codes  should 
xlso  be  modified  to  include  a  numerical  interpolation  scheme  to  generate  flexibil¬ 
ities  which  lie  between  the  tabulated  computational  design  data  points  contained 
herein.  Because  the  closed  form  advanced  strength  of  materials  beam  theory  is 
quite  accurate  for  many  cylindrical  geometries,  the  computer  codes  should  be 
modified  to  allow  the  code  to  choose  between  using  the  A.S.M.  beam  closed  form 
or  the  new  shell  flexibility  design  data.  With  these  design  data  shell  flexibilities 
included  in  the  transfer  matrix  analysis,  the  accuracy  and  range  of  application 
for  turbomachinery  is  greatly  improved. 
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APPENDIX  A 

Derivation  of  Equations  (1)  and  (2) 

Referring  to  the  general  shaft  element  of  Fig.  1,  the  load-displacement  rela¬ 
tionships  can  be  developed.  By  first  summing  forces  in  the  y  direction, 


^-i-v;  =  o  ,  v;  =  (.41) 

and  by  summing  moments  about  the  right  end, 

M{  =  Mi-i  Vi-il  (-42) 

By  solving  for  the  radial  displacement  at  the  right  end, 

Yi  =  V',-1  V-  9^1  -  Mt(fye)  -  Vdfyy)  (A3) 

and  for  the  angular  displacement  at  the  right  end, 

Si  =  0i- 1  ~  Mi(M  -  vdM  (A4) 

By  substituting  (Al)  and  (A2)  into  (A3)  and  (A4)  respectively, 

Yi  =  -  0i.il  4-  Mi-iiM  A  Vi-iKfy»)  ~  Yi-dfyy)  (A5) 

6,  =  9i_i  -  Mi-Moo)  *  V^lifee)  -  (fey)  (-46) 


where  fye  is  defined  as  the  flexibility  of  the  structure  for  the  radial  deflection 
under  bending  moment  loading.  fyy  is  defined  as  the  flexibility  of  the  structure 
for  the  radial  deflection  under  shear  loading,  fee  is  defined  as  the  flexibility  of 
the  structure  for  the  rotational  deflectior  under  bending  moment  loading.  fgy 
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is  defined  as  the  flexibility  of  the  structure  for  the  rotational  deflection  under 
shear  loading. 

Assembling  equations  (Al),  (A2),  (A5)  and  (A6)  into  matrix  form  to  be  used  in 
computational  analysis,  equation  (2)  of  the  main  text  is  derived, 


(2) 


Now,  by  defining  the  transfer  matrix  flexibility  coefficients  as; 

a  =  l(fy9 )  -  fyy 

b  =  fye 

<*  =  Kfso)  -  fey 
3  =  fee 

We  can  by  substitution  show  the  general  form  of  the  field  transfer  matrix  in 
terms  of  the  defined  flexibility  coefficients  which  is  equation  (1)  of  the  main 
text, 
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Assembly  of  Equation  (3) 

For  the  Timoshenko  beam  element,  the  flexibility  equations  derived  by  ad¬ 
vanced  strength  of  materials  are; 

l2 


fye  —  fey  — 
fyy 


2EI 

l3 


l 


fee  £I 


3  El  GACK, 

l 


(-48) 
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where  the  shear  form  factor  is 

6(l  +  i^)(l  ^  ^)2 

h‘  =  (7 T 6rv)( I  -  ^)2  —  {20°+  12i/)(^)2  ('49) 

By  substituting  the  expressions  of  (A8)  into  the  expressions  of  (A7)  and  sub¬ 
stituting  again  into  equation  (1),  the  general  field  transfer  matrix  for  a  beam 
element  becomes  equation  (3)  of  the  main  text, 


'1  i  JL  Jl _ l 

1  1  2  El  6  El  GACK, 

0  1^  -ii- 

U  1  El  2 El 

0  0  1  l 

.0  0  0  1 


Yi-t 

0,-1 

Mi-i 

Vi-i 
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APPENDIX  B 


Derivation  of  Equations  (4)  through  (10) 

Recalling  equation  (1)  for  the  general  shaft  element  of  Fig.  1, 

1 


y, 

9; 


v'  I 


■1 

l 

b 

a  n 

r  V'-! 

0 

1 

3 

a 

0 

0 

1 

l 

\ 

Mi-i 

.0 

0 

0 

1  . 

Vi-i 

(1) 


By  applying  the  cantilevered  boundary  conditions  used  for  multiple  station 
transfer  matrix  analysis; 

yt~i  =  o 

(B 1) 

Oi- i  =  0 

and  by  first  applying  a  shear  tip  load,  V',,  with  no  bending  moment  tip  load,  M, , 
the  force-displacement  relationships  can  be  written.  First  summing  the  forces 
in  the  y  direction, 

rt_i-Ft  =  0  ,  Vi  =  V) -i  ( 52 ) 

and  by  summing  moments  about  the  left  end, 

Mi-i  -VJ  =  0  .  = -\]l  ( 53 ) 


Now  substituting  (B2)  and  (B3)  into  (1)  and  solving  for  the  displacements 

V  j  =  aVt  -  bV,l 
9,  ~  a  Vj  —  3  Vt  l 

and  solving  for  the  flexibilities  a  and  a, 

y, 


(54) 


a  = - bl 

4-, 

*  i 


(55) 
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For  a  unit  shear  load. 


a  =  Y,  4-  bl 


a  =  9i  4-  31 

Bv  defining  1,,  and  0V  as  the  radial  and  rotational  displacements  respectively  of 
the  right  end  due  to  a  shear  tip  load  and  no  moment  load,  then  the  flexibilities 


become 


a  =  Yv  -r  bl 


a  =  9V  4-  31 


Now  applying  a  bending  moment  tip  load,  A/,,  with  no  shear  tip  load,  V, ,  the 
force-displacement  relationships  can  be  written.  First  summing  the  forces  in  the 


v  direction, 


V,-1  -1,-0  ,  Vj  _  i  —  V,  —  0 


and  by  summing  moments  about  the  left  end, 


Mi-i  -  M,  =  0  ,  Mi-i  =  Mi 


Now  substituting  (B8)  and  (B9)  into  (1)  and  solving  for  the  displacements, 


Y,  =  bM i 
9,  =  3.\It 


(B 10) 


and  by  defining  Ym  and  9m  as  the  radial  and  rotational  displacements  respec¬ 
tively  of  the  right  end  due  to  a  bending  moment  tip  load  and  no  shear  tip  load, 


Ym  =  bM, 

0m  =  3Mi 

For  a  unit  load,  equations  (5)  and  (7)  of  the  main  text  are  written  as 


(511) 


b  =  Y, 


3  =  0r 
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By  substituting  (5)  and  (7)  into  (B7),  equations  (4)  and  (6)  are  written  as 

(4) 
(6) 


a  =  l(Ym)  +  Yv 
a  =  l(0m)  +  0V 


(10) 
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APPENDIX  C 

Derivation  of  Equations  (11)  through  (14) 

Recall  the  general  structure  element  of  Fig.  1.  By  fixing  the  opposite  end  (i) 
and  loading  the  new  rigid  end  (i  —  1)  so  that  the  same  load  end  conditions  exist 
on  the  general  structure  element  as  before,  then  the  same  strain  energy  state 
exists  and  the  following  analysis  can  be  applied.  (See  Fig.  70  for  an  illustration 
of  the  general  structure  element  now  in  the  negative  attitude.) 

By  summing  forces  in  the  v  direction, 


1 —  V,  —  0 


Vt-x  =  1) 


and  by  summing  moments  about  the  right  end, 


(Cl) 


Mi-i  -  l(Vi)  -  Mi  =  0  ,  Mi-x  =  Mi  -  l(Vt)  (C 2) 


Solving  for  the  displacements  at  the  right  end  (t  -  1), 


Yi-i  =  Y  -  l(9t)  +  Mi.x(fy$)n  +  V-_a(/vy)n 
Bi-i  =0i-M,-1(Mn-Vl-1(fey)n 


(C3) 


where  “n”  denotes  the  flexibility  when  loaded  in  the  negative  attitude.  By 


assembling  into  matrix  form  and  isolating  the  flexibility  coefficients  in  order  to 
obtain  the  flexibility  transformation  matrix, 


Recalling  for  a  general  structural  element  that; 

a  =  l(  f  y$)  ~  fyy 


*=/,• 

ot  —  l{foe)  - 


(.47) 


■3  =  fee 
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substituting  “n”  for  the  negative  attitude  structure  but  loaded  as  in  Fig.  1,  and 
solving  for  the  flexibilities, 

{fy())n  —  in 

{fdd)n  —  3n 

(C5) 

( fyy)n  =  —  an 

( fdy)n  =  /(dn)  On 

By  substituting  (Co)  into  (C4), 


1  - 

0  1 


- 1  J  \ 
1  J  1  Bi  J 


Solving  for  the  displacements  at  end  ( i ), 


\9i  J 


in  l(i  n)  ®n 

~3n  Otn  —  l(^^3n'j 


f(i  n) 


A/i-i 

V'-! 


A/,-! 

Vi-i 


IM  [0  [0  lj[-dn  on  — /(/3n)J  |  Vi-,  J  '  ' 

Recalling  the  transformation  of  the  positive  attitude  general  structure  element, 
equation(l),  and  separating  the  flexibility  coefficients  in  order  to  obtain  the 
flexibility  transformation  matrix, 


i  / 1  r  rt_!  i  r  ip  aJ/Mi-O 

0  lj  \  9i-x  j  +  [3P  ap  J  \  Vi-1  J 


where  “p”  was  inserted  to  designate  the  flexibility  coefficients  of  the  positive 
attitude  general  structure  element.  Now  subtracting  (C7)  from  (C8), 


3p  otp 


1  / 1  -in  an  —  l(bn) 

0  lj[dn  Ki3n)-*n 


Solving  for  the  negative  attitude  flexibility  coefficients  in  terms  of  the  positive 
attitude  flexibility  coefficients, 

an  =  ap  -  /(ip)  -f  /2(/dp)  -  /(ap) 

i  n  /(dp)  ip 
On  —  /(dp)  Op 


3  n  —  3p 


(CIO) 
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These  equations  represent  the  transformation  of  the  flexibility  coefficients  of  the 
general  structure  element  from  the  positive  attitude  to  the  negative  attitude. 

Now  for  the  specific  case  of  the  conical  structure  used  in  turbomachinery, 
because  the  cone  is  an  axisvmmetric  structure,  the  flexibility  matrix  is  symmet¬ 
ric,  thus 

(Cll) 

then  by  substituting  (Cll)  into  (A7), 

a  =  13  —  b  (12) 

Also, 

OCn  =  l((3n)  - 

(cm 

ap  =  1{3P)  —  bp 

and  finally  by  substitution  of  (C13)  into  (CIO),  the  flexibility  coefficient  trans¬ 
formation  from  the  positive  attitude  of  a  conical  element  to  the  negative  attitude 
of  the  same  conical  element  is  written  as  equations  (11)  through  (14)  of  the  main 


Cln  —  &p 

(ii) 

bn  =  &  p 

(12) 

Qn  —  bp 

(13) 

3n  —  3p 

(14) 

text, 
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APPENDIX  D 


* 

ANSYS 

EXAMPLE  INPUT  FILE 

* 

EDWARD 

ANTHONY  L' ANTIGUA 

« 

TEXAS 

ASM  UNIVERSITY 

* 

* 

/PREP7 

/TITLE. ROTOR  WITH  SHELL  ELEMENTS 


C 

C  SELECT  SHELL  ELEMENT 

C 

ET, 1,63 
C 

C  SELECT  MATERIAL  PROPERTIES 

C  YOUNG'S  MOOULUS . POISSON ' S  RATIO 

C 

EX .  1  . 30E6 
NUXY ,1.0.3 
C 

C  INPUT  STRUCTURE  GEOMETR  ■' 

C  GENERATE  KEYPOINTS 

C  GENERATE  LINE  SEGMENTS 

C  GENERATE  AREAS 

C 

R.  1  ,3.2 

K, 1 ,0. 16 

K, 2. 4. 30. 9282 

K,  3, 8, 45. 8564 

L,  1  ,2 
L.2.3 

K, 1000,0.0 
K, 1001 .1,0 

AROTAT,  1,2 . 1000.  1001 .360,4 

C 

C  MESH  GENERATION 

C 

ELSIZ. .6 
AMESH. ALL 

C 

C  FIX  ALL  LEFT  ENO  NODES 

C 

NRSEL.X.O 
0 . ALL . ALL 
NALL 

C 

C  CREATE  INFINITELY  RIGID  RIGHT  END 

C  BY  USING  O.O.F,  CONSTRAINT  EQUATIONS 

C 

CERIG.  1  , 2 
CERIG.  1  .  3 
CERIG. 1 . 4 
CERIG.  1  . 5 
CERIG.  1  .6 
CERIG. 1 , 7 
CERIG. 1 . 92 
CERIG. 1 ,93 
CERIG. 1 ,94 
CERIG.  1  ,95 
CERIG. 1 .96 
CERIG. 1 , 97 
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CERIG. 1,170 
CERIG, 1.171 
CERIG. 1.172 
CERIG. 1,173 
CERIG. 1.174 
CERIG. 1.175 
CERIG. 1 . 248 
CERIG, 1 .249 
CERIG, 1 . 250 
CERIG, 1.251 
CERIG. 1 , 252 
C 

C  APPLY  RIGHT  END  LOADING 
C 

NRSEL.X.8 
NRSEL. Y.O 
F.ALL.FY.-0.5 
NALL 

C 

C  REORDER  SOLUTION  WAVEFRONT  IN  ORDER  TO 
C  MINIMIZE  MATRIX  SIZE  AND  SOLUTION  TIME 

C 

WSORT . X 

WSTART 

WAVES 

C 

C  ONE  ITERATION  FOR  LINEAR  SOLUTION 
C 

ITER. 1 
AFWRIT 
FINISH 
/EXE 

/INPUT. 27 

FINISH 

/EOF 

C 

C  POST  PROCESSING  DATA  COLLECTION 

C 

/POST  1 
/OUTPUT. 6 
SET .1.1.1 
C 

C  SELECT  MIDLENGTH  NODAL  DISPLACEMENTS 

C  TO  CHECK  FOR  ANY  OVALIZATION 

C 

NRSEL.X.4 
NRSEL , Y.O 
PRNODE 
PLDISP 
NALL 

4RSEL . X  .  4 

NRSEL. Z.O 

PRNODE 

PLOISP 

NALL 

C 

C  SELECT  END  NOOAL  DISPLACEMENTS  TO 

C  DETERMINE  FLEXIBILITY  COEFFICIENTS 

C 

NRSEL, X, 8 
NRSEL .Y.O 
PRNODE 
PLDISP 
NALL 

NRSEL  .  X  ,  8 
NRSEL. Z.O 
PRNOOE 
PLOISP 
NALL 

/OUTPUT.  1 
/EOF 
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